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STABILITY OF PLANE HA GNETOHYDROPYNAMIC CHANNEL FLOV7 
. I'JITH PARALLEL MAGNETIC FIELD 

P. R« Nachtsheim,’* MASA-Lewis Research Center 
end 

E, Reshotko, - Case Institute of Technology 
Cleveland, Ohio 

An approximate solution to the title problem was given in a paper by 

J. T. Stuart‘S where a reduced fourth-order disturbance .differential equation 

proposed for small -magnetic Reynolds number was solved using asymptotic 

) 

methods. The problem has been reexamined through exact numerical integra- 
tion of the pertinent sixth-order system of disturbance equations throughout. 

2 

an extended range of magnetic Reynolds numbers . For comparison, exact 

numerical results were .also obtained to Stuart's fourth-order equation. The- 

results indicate that Stuart's simplifying assumption is justified only for 

magnetic Reynolds numbers small compared to one. For magnetic Reynolds 

numbers of order one- or greater, there are significant changes in the sta- 

% 

bility characteristics reflecting the- increased ‘importance of magnetic 
effects. This is also borne out in a calculation of the various viscous and 
magnetic contributions to the rate of change of disturbance energy. It is 
shown that the resistivity enters this problem in two ways: (l) it sets up a 
time-independent Maxwell stress that augments the disturbance energy when it 
is 'of the same sign as the vorticity of the basic flow, and (2), through joule 
dissipation the disturbance energy is decreased. For small resistivity 
(Re^» 1) the former dominates, leading to a net augmentation of the disturbanc 
energy, while for large resistivity (Re 1 ) the dissipative effect is 
dominant. 

Now at NASA-Ames Research Center 

1 Stuart, J.T.- : Proc. Roy Soc A 2 21 18? (1 9$h) 

2 

Nachtsheim, P.R. and Reshotko, E,: NASA TN December, 
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MAGNETOHYDRODYNAMIC BOUNDARY BAYERS 
INVOLVING NON-EQUILIBRIUM IONIZATION 


A, Sherman 

General Electric Space Sciences Laboratory, Valley Forge, Pa. 

and 


E. Reshotko 

Case Institute of Technology, Division of Engineering, Cleveland, O. 
(Consultant, General Electric Space Sciences Laboratory) 


In the study of practical magnet ohvdro dynamic devices such as 
generators and accelerators recent experimental and theoretical 
research has shown the importance of non -equilibrium ionization for 
their successful operation. At the same time, these experimental 
results have demonstrated the existence of large voltage drops at the 
electrodes', thereby preventing full utilization of the non -equilibrium 
ionization. These voltage drops- occur not only in the thin sheath at the 
base of the channel- boundary layer, but also partially in the bo un dary 
layer itself. Presently available theoretical analyses of the magneto- 
hydrodynamic boundary layer take into account non- equilibrium 
ionization in' only an approximate manner, and neglect the sheath as 
well. They have not been able to predict the large voltage drops 
observed experimentally. 

The present paper will formulate the magnetohydrodynamic 
boundary layer problem' when non -equilibrium ionization exists. In 
particular it will be shown that when non -equilibrium ionization is taken 
into account properly, that the characteristics of the sheath at the base 
of the continuum boundary layer are essential to the correct solution of 
the problem. In addition, since non -equilibrium ionization occurs more 
readily when segmented electrodes prevent Hall currents, it is Important 
to be able to analyze the boundary layer over an electrode wall with 
alternate finite electrode and insulator segments. For this reason 
solution- of the equations will be carried out by a finite difference technique. 



PROBLEMS' H JxB PLASMA ACCELERATION 
James A. Fay, Marvin Goldstein and Peter Sockol 
Department of Mechanical Engineering 
Massachusetts Institute of Technology 

This research is concerned with problems encountered in accelerating 
collision-dominated plasmas by an externally applied magnetic field inter- 
acting with currents also supplied from an external source. A principle 
focus of the research is the interaction of the current-carrying plasma 
with solid boundaries* Summarized briefly below are accounts of progress 
on four related problems. 

(a) Experimental measurement of heat transfer to the walls of a 
Faraday shock tunnel is being investigated. Heat transfer is measured by 
monitoring the temperature rise of "thin" and "thick" film heat transfer 
gages with an infrared detector. Two types of gages are being developed. 

The first is a heat capacity gage with 0.1 microsecond response time for 
measuring rapidly varying heat transfer. The gage consists of carbon and 
aluminum coatings applied to thallium bromide or arsenic trisulfide windows, 
the temperature being monitored by a mercury doped germanium detector at 
liquid helium temperature. The second gage which is being developed to 
measure heat transfer to an emitting electrode , consists of a thicker 
metallic coating for which the response time is about 10 microseconds but 
which is sufficiently thick to conduct the emitted current without being 
heated by Joule losses in the electrode surface of which the gage is an 
integral part. At the present time problems of gage coating, infrared • 
transmission, calibration, etc. are being studied. Also, a test section for 
use with the gages is being tested. 

(b) A new formulation for transport effects in multi-component 
boundary layers has been devised which shows promise of simplifying the 
numerical solution of plasma boundary layer problems. This new method in- 
volves .direct use of certain moments of the linearized Boltzmann equation 
from which the gradients of conserved properties are linearly related to the 
fluxes of these same quantities, and bears some resemblance to the transport 
limit of the Grad moment method. In combination with the conservation laws 
relating the divergence of the fluxes with convective derivatives and source 
terms, a set of first order differential equations is obtained (total dif- 
ferentials in the case of similarity-type boundary layers) In which the usual 
multi -component transport coefficients , such as the total thermal conductivity 
do not appear explicitly. An additional advantage of this method is that 
different temperatures of electrons and heavy species and different levels of 
approximation in the Sonine polynomial expansion of the perturbation to the 
distribution function may be used for different species, which is important 
for plasma boundary layers. A report is being written for distribution. 

(c) A theoretical study of Hartmann boundary layers in a homo- 
polar device has yielded preliminary results which are in approximate 
agreement with experimental measurements. The theory takes into account two 
energy equations, (one for the electrons and a second for the heavy particles) 
a diffusion equation for the electron-ion pairs and a momentum equation for 
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the heavy species. Approximate solutions, which use average values of the 
transport coefficients, show that the flow velocity in the homopolar 
device, is an appreciable fraction of the "ionization- velocity," i.e. , the 
velocity of an atom whose- kinetic energy equals its ionization energy, and 
that appreciable variation in the. degree of ionization can result from 
.small changes in the total current, at least for the small degrees. of ioniza- 
tion for which the present approximate solutions are valid, 

(a) A zero-order theory for the stabilization by a magnetic field 
of an arc transverse to a flowing steam has- been devised. It is based on_ a 
model in which the hot arc column is motionless, but is in static equilibrium 
with the dynamic pressure of the "cold" stream in which it is immersed, and 
to which heat is conducted from the hot column. In this zero-order theory, 
the cross-sectional shape of the arc column is determined from the momentum 
equation and the temperature distribution within the column is determined 
from the energy equation. The solution of these equations leads to relation- 
ships between the stabilizing magnetic field strength and the flow- velocity, 
between the electric field and the width of the column, and. also between the 
arc power, and the average temperature. The first twg relationships are in 
approximate agreement with the measurements of Myers . 


Meyers, T. , "Experimental Investigation of the Balancing Mechanism and 
Electrical Properties of an Argon Arc in Crossed Convective and Magnetic 
Fields," Tth Symposium on Engineering Aspects of MKD. (19 66) 
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ELECTRODE EFFECTS IN ACCELERATORS 

By 

S . Aisenberg 
P. Hu 
V. Rohatgi 
S. Ziering 

SPACE SCIENCES, INC. 

301 Bear Hill Road 
Waltham, Massachusetts 


1 . Experiments and Analysis 

A number of problems are investigated experimentally and analy- 
tically in order to give a better understanding of some of the important 
physical processes involved in plasma surface interactions. The 
features studied include: the tangential drag on electrodes in J x B 
accelerations; the arc construction processes; and a theory proposed 
to explain electron emission from cold cathodes. by an ion microfield 
emission process. 

The tangential electrode drag forces are measured for an arc in a 
transverse magnetic field, and it is shown that the cathode and the 
anode tangential forces are an appreciable fraction of the J x B force 
driving the plasma. The presence of tangential electrode (and boundary) 
friction forces is expected for a hot plasma flowing in a channel when 
one considers that some of the tangential momentum is transferred each 
time an ion, electron, or gas atom collides with the surface. 

The basic features of the electrode structure used in this experiment 
are shown in Figure 1. The force measurements are made for the central 
electrode which is operated as the cathode or alternatively as the anode. 
This central electrode is supported on a torsion spring structure so that 
the tangential force (or torque) can be determined from the angle of 
rotation . 

Some results are illustrated in Figure 2 where the measured tangential 
force is shown as a function of magnetic field for a 77 amp dc arc in 
24 Torr of argon, and where the central electrode is operated as a cathode 
as well as an anode. The cathode force due to ion current drag is 
deduced by assuming that the neutral gas component at the cathode is 
the same as the anode drag. (This assumption should be examined 
further because the boundary layers at the anode and cathode may be 
very different) . The measurements show that: 
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1) There is a pronounced threshold in the cathode and anode force as 
a function of magnetic field and of arc current. 

2) If it can be assumed that the momentum transfer to the anode is 
predominately due to the neutral gas atoms , there appears to be a 
threshold in the electrode drag by the neutral gas components . 

3} The cathode drag appears to show a saturation. 

4) The cathode force is as large as 24% of the input force driving the 
plasma, and the anode force is. as large as 13%. 

5) ' The deduced ion current component is larger than the neutral gas 

component, and in the range measured, is about 16%' of the plasma 
driving force. 

The tangential force on an electrode is proportional to the tangential 
velocity, to the perpendicular particle current flow , and to the tangential ' 
momentum accommodation coefficient. The observation that the tangential 
force at the cathode was larger than the anode force indicates that the’ 
positive ion current flow to the cathode is capable of transferring an 
appreciable tangential force to the, cathode. The direct measurements of 
the electrode drag forces show that the effect is serious enough to be 
considered further and at the same time gives- interesting data that 
should be considered by any proposed models . 

An analysis was made of, the arc constriction processes because the 
constriction influences- the plasma surface interaction. The constriction 
mechanisms considered involve the balance between the input power to 
the plasma and the following modes of power loss: 1) diffusion loss for, 
low pressure plasmas, 2) radiation loss for high density, high temperature 
plasmas, 3)- thermal energy loss for higher pressure plasmas, 4) convec-r 
tive energy loss due to the flow of cooler gas. The available data and 
the theory both show that for the diffusion and for the thermal conductivity 
mode loss, the product of the current density T and the arc current 
should be constant. 

2 . Theory 

In order to clarify plasma electrode interactions and 'facilitate a 
clearer understanding of the physical parameters governing sheath 
phenomena, the existing collisionless theory of a plasma sheath near . 
an infinite planar electrode is examined in detail. Ranges' of validity, 
as well as difficulties with the existing sheath models are investigated. 
The assumption of monoenergetic particle distributions, as used for 
instance in “Bohm's Sheath Criterion” restricts the usual analysis. To 
remove this difficulty Maxwellian particle distribution functions are 
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assumed. In addition, microscopic boundary conditions are assumed 
to incorporate the respective electrode and plasma parameters . A new 
sheath model is then developed without imposing the usual sheath 
criterion for the conditions specified at the sheath edge. Explicit 
results for the current-voltage relation and for the sheath thickness 
are obtained for the improved sheath model proposed in this study. 
Other important physical quantities are expressed in terms of the 
electrostatic potential, which can be calculated by a quadrature sub- 
ject to the physical parameters to be specified. The general solution 
reduces to a simple analytic form for the special case where the tem- 
perature of the electrons , ions and the electrode are identical. 

Based on the above work, the next logical step in a systematic 
analysis necessitates consideration of ionization processes in the 
domain of the pre-sheath as well as the sheath itself. Towards this 
end it was found necessary to construct kinetic equations which in- 
corporate ionization processes. Ionization processes are usually 
introduced into the equations of magnetohydrodynamics as gross para- 
meters to be determined by experiment. The mathematical accessibility 
of model kinetic equations in neutral and ionized gases, motivates the 
formulation of kinetic equations for a three-component plasma which 
incorporate the ionization process in this study. As a first attempt, 
only the ionization processes resulting from electron-neutral collisions 
are considered. The resulting mathematical formalism yields a set of 
model equations which are mathematically not more complicated than 
those for plasmas without consideration of ionization. The possible 
application of the present model equations to various problems is 
discussed, such as the ionization growth in a gas, the cathode sheath, 
and the structure of an ionizing shock wave . 

Reports and Publications 

The following publications have been generated as a result of the 
research performed under this contract; 

1) Aisenberg, S. and Rohatgi, V; 

"Measured Tangential Electrode Forces for an Arc in a Transverse 

Magnetic Field". App. Phy. Letters, 8, 194 (1966) 

2) Aisenberg, S. and Rohatgi, V.; 

"A Study of Arc Constriction Processes", Presented at the 7th 

Symposium on Engineering Aspects of Magnetohydrodynamics in 

Princeton, N.J.. , March, 1966 
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3) 'Aisenberg, S. and Rohatgi , V*; 

"A Metallic Insulator for Arc Electrodes' 1 , (to be submitted for 
publication) 

4) Hu, P.N., and Ziering, S., "Gollisionless Theory of a Plasma 
Sheath Near an Electrode " (submitted for publication) 

5) Hu, P.N. and Ziering, S., "Kinetic Model for Three ^Component 
Plasmas With Ionization" 1 , (submitted for publication) 
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MAGNETOAERODYNAMI C DRAG AND 
SHOCK STAND-OFF DISTANCE 


C. F. Chang, W R. J. Nowak, ^ S.' Kranc, ^ R. Porter ^ 

G. Trezek, ^ M, C. Yuen, ^ T. P, Anderson, ^ and A. B. Cambel^ 


Gas Dynamics Laboratory 
Northwestern University 


It has been shown that the incorporation of magnetic field coils in * 
space vehicles offers interesting, possibilities in alleviating some of the 
problems of space flight. For example, the coupling of the flow field with 
a magnetic field has the salutory effect of increasing the drag experienced 
by an entry vehicle. Other promising applications of external magnetogas- 
dynamics are: the opening of communications windows, enhanced flexibility 
of maneuvers by virtue of flight control, active shielding against radiation, 
the reduction of convective heat transfer and the elimination of the non- 
reusable heat shield. So far, such schemes have not been exploited due to 
hardware design limitations and due to the fact that current design philoso- 
phy has been successful in current programs. However, continuing progress 
in superconducting magnets and the demands of advanced space programs lend 
practical significance to equipping space vehicles with powerful electro- 
magnets. The purpose of this program is to investigate this problem 
systematically from both analytical and experimental viewpoints. . 

A theoretical study was initiated in order to study and compare the 
factors important in blunt body, flight magnetohydrodynamics and laboratory 
simulation in the thermal arc plasma facility. 

-Consideration of the flow regimes indicates that assumptions of 
continuum flow, a Hugoniot shock, and sub- layers of viscous and sheath effects' 
are valid for entry flight below 250 kilofeet. The assumptions of a Hugoniot 
shock and thin boundary layer are considerably less valid for the laboratory 
conditions owing, to the low Reynolds numbers. Therefore, a viscous shock 
layer was included in some of the analyses discussed below. Theoretical 


(1) Graduate Students 

(2) Faculty 

(3) Chairman and Walter P. Murphy Professor of Mechanical Engineering and 
Astronautical Sciences. 
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difficulties prevented including the shock structure, however. 

The effects of viscosity, deformation of the magnetic field and the 
Hall effect on the magnetic interaction were studied by the local similarity 
method for the stagnation region of aerodynamic - like flow. The flow 
observables were found for a range of dimensionless parameters appropriate 
for the flight and laboratory conditions. 

It was found that the effect of low Reynolds number (laboratory flow) 
is very important for low magnetic interaction parameter but that a boundary 
layer develops at high interaction# Thus, the outer flow becomes essentially 
inviscid for large Hartmann number similar to MHD Poiseuille flow. The 
effect of the deformation of the magnetic Reynolds number was negligible for 
magnetic Reynolds numbers up to an order of magnitude greater than that for 
unseeded, equilibrium plasma. However, an increase of an additional order 
of magnitude resulted in the field being considerably distorted with a 
wrapping of flux lines about the body. The stand-off was decreased and the 
drag increased for increasing magnetic Reynolds number. The analysis showed 
that it is incorrect to assume that a dipole source looks like a dipole at 
the body surface when there is distortion but that the assumption of a 
(reduced) dipole strength at the shock is valid. The Hall effect is very 
important for flight and laboratory conditions with kilogauss field .strengths 
at the stagnation point. The magnetic interaction is reduced more so for 
an insulated body than a conducting one but the effect is important for 
both. A rolling moment is introduced because of azimuthal flow generated 
by the Hall effect. All of the special effects tend to increase the critical 
interaction parameter associated with magnetic support where the attached 
shock layer theory fails. 

Experimental investigations of magnetoaerodynamic drag have shown 
substantial (~15%) increase in drag with magnetic field. Tests have been 
run under a variety of flow conditions with both 1-1/2” and 3" bodies. In 
general the experimental aerodynamic drag component agrees reasonably with 
theoretical predictions. The measurements of drag with field for the 3” 
body do not agree with prediction nearly so well as do those tests for the 
1-1/2” body. The main reason is the non- uniformity in free stream conditions 
which affect the 3” body more than the 1-1/2" body. 

Experimental values of drag increase are compared to theory for a 
regime where the viscous and Hall effects are important. No theory exists 
for the combined effects. The magnetic interaction parameter is obtained, 
from experimentally measured quantities except for the electrical conducti- 
vity which is estimated from equilibrium calculations. We have found 
that while drag appears to be greater under conditions where the Hall and 
viscous effects are lessened, the various data may be brought closer by 
plotting against the product of the interaction parameter and the shock 
density ratio. The remaining differences may be attributed to the Hall 
effect. Future work will include the measurements of different drag 
components. 
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The shock stand-off distance is plotted versus magnetic field squared. 
These measurements were made by analyzing photographs of the flow using a 
scanning microdensitometer . We have found that the experimental measure- 
ment is severely hampered by a thickened shock under our running conditions. 
Nevertheless considerable alteration in the shock stand-off has been observed. 
Present efforts include the measurements of electron density and temperature 
across the shock structure spectroscopically in order to better define the 
shock stand-off distance. 

In addition to the experimental drag studies, a program has been 
initiated to experimentally determine magneto aero dynamic heat transfer to 
a blunt body. A thin shell, 3' 1 diameter, hemispherical, copper model has 
been instrumented with thermocouples . In the presence of an argon plasma, 
the local net transient heat transfer rate to the model is measured with, and 
without the presence of a magnetic field. For a thin shell the heat 
transfer rate to the model is proportional to the temperature - time rate 
(dT/dt) of the hody if temperature gradients in the body are neglected. 
Preliminary results with a magnetic field strength of 7700 gauss at the 
stagnation point yield a 26% reduction in model centerline dT/dt. However, 
corrections must be made to account for temperature gradients which do 
exist in the body during the test time. 

Future tests will make possible the mapping of the magnetoaerodynamic 
heat transfer over the surface of the body (to half angles of 45°) for a 
range of magnetic field strengths of 8000 gauss and lower. 



PUBLICATIONS 
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Properties of a Magnetically Suspended Arc in Supersonic Flow 

by 

Arnold M. Kuethe, University of Michigan 

A NASA grant on which work began March 1, 1966, covers experi- 
mental and theoretical studies of the properties and structure of a dc arc 
suspended in mutually perpendicular electric, magnetic, and airflow 
fields. Previous work under an Air Force contract, in collaboration 
with Charles E. Bond, showed 1) that the arc can be stabilized on rail 
electrodes at a slant angle near the Mach angle of the flow, 2} that the 
'•aerodynamic drag” of the arc changes with arc current in such a way 
that the magnetic field required to stabilize the arc is constant over a 
wide range of the arc current, 3) column rather than root phenomena 
determine the conditions for arc stability, and other gross properties 
of the arc. 

The work under the NASA grant is being directed toward the 
details of the fluid mechanical structure of the arc and the transfer 
processes by which mass, momentum, and energy are transferred 
to the ambient flow. 

Measurements to date indicate that, over the range of the 
measurements {Mach numbers up to 3. 5, ambient pressure down to 
10 mm Hg) the arc simulates a solid cylinder with a drag coefficient 
between 1. 0 and 1. 5 balancing the Barents force. The solid body 
hypothesis is based primarily on photographic observations of an axial 
flow from anode to cathode of about 10 meters/sec in the Mach 2. 5 arc. 

In addition to the axial flow, the incremental Borentz force on any 
relatively hot filament, (i. e. , higher conductivity and current density 
than the surroundings) will force it upstream to the arc boundary where 
it will be cooled and carried downstream with the flow along the boundary. 
This circulation will resemble that caused by a vortex pair, providing 
an upstream flow toward the stagnation point. Thus a model of the arc 
in agreement with the measurements and physical reasoning is as follows: 
The streamlines projected on a cross-section would resemble those of a 
distributed source at the center, a vortex pair and a distributed sink at 
the boundary. 

The work will cover the widest possible range of Mach numbers, 
ambient pressures, arc current and magnetic field. Wake traverses 
of gas properties will be, made. Theoretical analyses will include arc 
stability and transport phenomena. 



(o) View looking downstream 



(b) View looking West 


Figure ?. Experimental Setup 
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A CRITICAL MASS FLOW MODEL FOR THE MFD ARCJET 


,H. A. Hass an 

North Carolina State University 
Raleigh, North Carolina 


.NsG-363 


At the recent Fifth Electric. Propulsion Conference a phenomenological 
theory describing the MPD arc jet was presented By Bennett, John, Enos and 
Tuchman. The theory employs the concept of an effective mass flow rate 2 >3 
and is Based on the assumption that the -arc discharge will operate at an 
-effective mass flow rate, m e , which minimizes the arc voltage independent 
of the actual mass flow rate passing through the device. As a result o£- 
this theory it is shown that, for mass flow rates greater or equal to 
the maximum, attainable velocity is limited by a critical velocity, v , 
given by v c - (21/m)^-' 2 where I is the ionization potential and in ie C the 
mass of the atom. For mass flow rates less than * e the maximum velocity 
can exceed v c . .It appears that the concept of effective mass flow rate, 
was introduced^ to explain- the existence of the voltage plateau which is- 
characteristic of rotating plasmas. Since it has been shown'*'’’ 5 that an 
interpretation of this phenomenon can Be obtained from a consideration of 
the conservation equations., the observed behavior of the MPD arc jet should 
be explained without invoking any assumptions' regarding the mass flow rate 
of the voltage. 

It has been established experimentally and theoretically that the volt- 
age" has the representation V = V 0 f CB where' V 0 and C are constants depending ‘ 
on the propellant- and B is the magnetic field strength. Using this relation 
and the conservation equations, it is shown that there exists a critical 
mass flow rate at which the thrust is maximum. In ;the absence of losses, 
the efficiency at this critical mass- flow rate is 50% and, for mass flow 
rates less (or greater) than the critical mass flow rate, the' efficiency 
is greater (or less) than 50%. The critical mass flow increases with an 
increase in total current and/or magnetic field strenth. This critical 
mass, flow rate plays the same role played by of Ref. 1. 

An experimental verification of the above theory requires thrust measure- 
ments for a given total current and magnetic field strength for a range of 
mass flow rates at low tank pressure. Such measurements were reported by 
Broclonan, Hess, Bowen and Jarrett 6 for various tank pressures. The figures 
show a least-square fit of the data of Ref. 6 for low tank pressures. These 
- figures indicate the existence of an optimum mass flow rate which increases 
with an increase in the total current and/ or magnetic field strength in agre 
meat with the predictions of the theory. 



References 


Bennett, S. , et al. "Experimental investigation of the MPD arc jet. " 

AIAA Paper 66-239 (March, 1966). 

^Moore, R. A. , Cann, G. L. , and Gallagher, L. R. "High specific impulse 
thermal arc jet thruster technology." AFAPL-TR-65-48, Part 1. Air Force 
Aero Propulsion Laboratory (June, 1965). 

3 

Cann, G. L. , et al. "Hall current accelerator," NASA CR-5470 (1966). 

4 

Hassan, H. A. , et al . "Experiments and analysis for coaxial Hall current 
accelerators and the role of ionization effects, " Sixth Symposium on 
Engineering Aspects of MHD, 83-84 (April, 1965). 

5 

Hassan, H, A,, Hess, R. V., and Grossmann, W. "Study of coaxial Hall 
current accelerators at moderate pressures." To appear as a NASA TN. 

6 

Brockman, P. , et _al, "Diagnostic studies in a Hall accelerator at low 
exhaust pressure." AIAA Paper 65-297 (July, 1965). 




g rams 



24 


ANALYSIS OF INSTABILITIES IN A LINEAR HALL CURRENT ACCELERATOR 


G. W. Garrison, Jr. and H. A. Hassan 
North Carolina State University 
Raleigh, North Carolina 

NsG - 363 


Experimental investigation of instabilities in linear Hall current - 
accelerators with applied axial, electric fields and radial magnetic fields . 
showed the existence of screw type instabilities^ >2 which show a striking 
resemblance to the instabilities observed in the positive column. Various 
attemptsl>2 have been made to apply the various - theories developed by 
Kadomtsev and Nedospasov^, Simorr*' , Hoh^ and Morse^ for the study of in- 
stabilities in the positive column and crossed field geometries to the 
linear- Hall current accelerator. However, none of the above theories was 
directly applicable. 

Because of the observed instability in the linear Hall current accelerator 
is similar to that observed in the positive column, the governing equations 
employed here are those normally employed in the study of weekly ionized gases. 
These equations were also the basic equations employed in References 3-6. A 
normal mode analysis similar to that employed by Johnson and Jerde^ for the 
positive column has been employed. It is shown that, regardless of the- direction 
of the., density gradient, no screw type instability exist if it is assumed that 
the axial magnetic field." is identically zero. Since the driving mechanicm .of 
the, observed instability is an E x B drift the above result is expected because 
a RjT component alone does not yield a radial drift component. 

Since the axial component of the magnetic field is small but never identi- 
cally zero in a linear 'Hall current accelerator, the analysis has been carried 
out for a magnetic field having both a radial and axial components. An exact 
solution of the governing equations resulted in a density distribution with no 
axial or azimuthal dependence. This solution was expressed in terms of the 
confluent hypergeometric function. Assuming that the perturbed density and 
potential have the general representation f(r) exp i ( u)t + "kg + m0) whereto 
is the frequency and It .is the wave number a general dispersion relation has been 
derived. The special case of a right handed helix (m - -1 mode) has been dis- 
cussed in detail and the results are presented graphically. The Figures show- 
plots of frequency, wave length, electric field for typical mobilities vs. the 
critical value of Bg/Bj. at the stability boundary. It is seen that the dimen- 
sionless parameter x c which depends on the ionization rate, radius of the device 
and radial and axial magnetic field strengths has a great influence on the stabil- 
ity! as the value of x c increase .the discharge becomes more and more unstable. 

As .an illustration, the calculations show that the discharge is unstable for 
values of x Q between 1,5 and 2.0 when 6 , the mobility ratio, is 10. This 

result Is somewhat interesting because the results of Ref. 1 indicate an onset 
of instability as soon as the magnetic field is switched on. The results also 
show that, at the stability boundary, the wave length is almost independent 
of the magnetic field strength except for a combination of high x Q and B K /B r . 

Also, an increase in Bj. promotes stability. 
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EXPEKTMENTAL INVESTIGATIONS OF STRONG SHOCK WAVES 
MOVING THROUGH AN IONIZED GAS 
G. L. Spencer 

Case Institute of Technology 
Cleveland, Ohio 


An experimental investigation of shock waves moving through a highly v 
ionized gas is reported. The shock waves are generated in an electromagnetic 
T-tube. The ambient gas, helium, is highly preionized by means of an axial 
current discharge with the shock waves being initiated during the after-glow 
period of the preionization. Measurements have shown that the plasma is in ‘ 
local thermal equilibrium and nearly completely ionized during the early 
after-glow period of the preionization. Typical temperatures and electron 
densities are in the order of 22,000°K and 3 x 10 16 electrons/cm 3 . Gasdyn- . 
amic shock waves were studied and the results compared to the theory presented 
by Jaffrin and Prob stein.* It was possible to verify the existence of, and 
quantitatively measure the thickness of the first layer, as predicted by 
their theory. Magnetohydrodynamic effects were studied via provision of an 
axial magnetic field. Magnetic field jumps and currents were observed in the 
shock front and found to be completely stable as long as the shock wave 
could be observed. It was found that these effects did not appear if either 
the applied m agnetic field or preionization were absent. 


*Jaffrin, M. Y. and R. F. Probstein, Phys. Fluids 7, 1958 (196*0. 



30 


HARMONIC GENERATION IN A MICROWAVE PLASMA 
F. J. Mayer and 0. K. Mawardi 
Case Institute of Technology 
Cleveland-, Ohio 


A number of attenpts have been made to estimate the harmonies generated 
in a microwave discharge.! These predictions have been compared with exper- 
iments . ^ The corroboration of the theory with the experimental observations 
has often been artifically improved by the arbitrary choice of parameters 
appearing in the theory. One such parameter is the electron temperature . 

Now as this temperature has a marked influence on the amplitude of the har- 
monics one would expect that the temperature should be found as one of the 
consequences of a correct theory. 

Accordingly, the affect of a few parameters were studied systematically. 
It was demons trated3 that the temperature dependence of the collision fre- 
quency was unimportant in as far as explaining the result of Uenohara ' s 
experiment. 

Work is now proceeding to investigate the effect of density gradients 
and ionization on the harmonic generation. 


X B. Murphy, Phys. Fluid 8, 1534 (1965); P. Rosen, Phys. Fluid A, 3^1 (1961). 
2 M. Uenohara, et.al. , Proc. Ionization Phenomena in Gases, 2 _> 768 (I960)., 

%. Mayer, Phys. Fluid (In Press). 
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A WIDE-BAND DICKE TYPE RADIOMETER 
A. T. Alper and 0. K. Mawardi 
Case Institute of Technology 
Cleveland, Ohio 


A wide-hand Dicke type radiometer capable of operating within 2.6 to 
14 KMC frequency range has been constructed. In order to achieve this ex- 
tremely wide bandwidth operation double ridge waveguide instruments, such 
as a balanced mixer, a directional coupler, a ridged-waveguide switch, two 
coaxial to ridge adapters, and X and S-band to ridge waveguide transitions, 
have been designed and developed. Difficulties have been encountered in 
achieving uniform frequency response of the balanced mixer due to matching 
of mixer T, E n and ’’H 1 ’ planes to the balanced arms. Aside from 4.5 - 5 KMC 
frequency bandwidth, the radiometer has been calibrated in the range of 
2.6 to 9 KMC using argon and fluorescent noise sources. The radiometer sen- 
sitivity, equivalent to a minimum detectable temperature change of 1.6°K 
has been computed. 

With wide-band capability of this radiometer, electron cyclotron radi- 
ation intensity measurements in a plasma have also been made by sweeping 
local oscillator of the radiometer from 2.6 to 9 KMC. 

The radiometer is planned to be used in a series of investigations 
dealing with electron cyclotron harmonic generation in a plasma. 
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EXPERIMENTAL ST®! OP THE INTERACTION OP ENERGETIC ELECTRON 

AND A PLASMA 

W. B. Johnson and M. R. £taith~ 

Case Institute of Technology 
Cleveland,' Ohio 


An experimental^ determination is made of the distribution of energy 
losses suffered by a" beam of energetic test electrons traversing a dense, 
hi#! temperature plasma. The test electron energy is 3,000 electron volts 
and the plasma density varies from 5 x 10 1 3 to 1 x 10 16 cm" 3 at a temper- 
ature of 'v 4b,000°K. The experiment was designed to produce data on single 
electron interactions .and not on beam-plasma effects. The experimental 
average energy loss is observed to be larger than the theoretical prediction, 
indicating a minimum discrepancy of at least an order of magnitude between 
theory and experiment. The theory includes the effects of binary collisions 
and collective interactions'. The discrepancy probably is explained by the 
neglect in the theory of higher order terms in the plasma kinetic equations. 
-The experimental average energy loss and spread in energy are both- found to 
vary, proportional to the plasma electron density. 1 * 2 


Reported in detail in: 

i 

X 

,! A Measurement of. Energetic Test Electron Interaction with a Plasma” , M. R. 
Smith and W. B. Johnson, Plasma Research Program T.R. A-38, June, 1965; 
Case Institute of Technology. 

^"Collisional Interactions of Energetic Test Electrons with a Plasma", M. R. 
Smith, Plasma Research Program T.R. A-39, July, 1965; Case Institute of 
Tprhnol rsirv; 
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A BEAT FREQUENCY INTERFERQMETER FOR PLASMA DIAGNOSTICS* 
W. B. Johnson, A. B. Larsen, and T. P. Sosnowski ' 
Case Institute of ' Technology 
Cleveland, Ohio 


A laser method has been developed for the diagnostics of plasmas with 
electron densities greater than about 10 11 electrons/cm 3 . A spatial resol- 
ution of about one microsecond. Is calculated for the particular apparatus 
which has been used. The range of sensitivity corresponds favorably to that 
achieved by microwave methods, but the spatial resolution is much superior 
to that afforded by microwaves. The method depends on the detection of the 
shift in frequency of a laser when a plasma Is introduced into the optical 
cavity. 

If a cavity of length L is employed to probe a plasma of length SL in- 
ternal to the cavity, then a frequency shift Av is expected which is given 
by 



CO 2 k 



Cl) 


where ur 2 = 4ime 2 /m is the square of the plasma frequency, v is the oscil- 
latory frequency of the laser, and n 0 is the index of refraction of all other 
particles in the plasma. To separate the effects of plasma, the frequency 
shift must be measured at two laser frequencies. Conveniently He-Ne lasers 
operate well at both 6328 A and 11,523 A. At the first wavelength, the ex- 
pected plasma contribution for a density of 10 11 electrons/cm 3 is 4.3 kc/s 
for 1/L = 1/2. Observed frequency shifts are about two orders of magnitude - 
higher, which shows the necessity for making the dual measurement. PM de- 
tector techniques have been used with frequency differences up to about 
5 mc/s. For larger shifts, spectrum analyzer methods could be employed out 
to the limit of photomultiplier tubes. Thus these techniques complement 
both microwave and other laser methods for the measurement of plasma electron 
densities. 


Experimentally a system has been constructed which employs two -laser 
cavities; one laser serves as the probe and the other the "local oscillator" 
laser. Each of the lasers is isolated optically and electronically but 
coupled together mechanically to assist in the cancellation of mechanical 
and thermal drifts. The output beams from each laser are made to coalesce 
by a beam splitter and are allowed to illuminate the surface of a photo- 
multiplier. Two channels are thus available to make measurements simultan- 
eously at the two wavelengths, 6328 A and 11,523 A. Using these techniques 
and others to be described in the paper, excellent overall system stability 
has been achieved in standard laboratory surroundings. Beat frequencies of 
* 150 kc/s for at least three minutes. The former corresponds to a stability 
of the total optical path length to about ± 0.1 A. Short-term stability 
measurements therefore show that plasmas in the order of 10 11 electrons/cm 3 
can be probed. Experimental data on argon discharge tubes will be presented. 



It is felt that with a better environment that it should be possible to im- 
prove on this stability by an order of magnitude . 


*Also published in: 

Journal of Quantum Electronics, April, 1966, Vol. QE-2, Paper 8C-8 Pg. lx. 
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TURBULENCE IN A RAREFIED PLASMA 
C. M. Tchen 

I 

, National Bureau of' Standards, Washington, D. C. 

A theory is developed on the damping process' in a turbulent 
plasma (nonlinear Landau , damping) . The Vlasov-Poisson* equations., 
describing the turbulent fluctuations, generate a chain of equations 
for their correlations. If the fourth order correlation. is degenerated 
into lower orders, the chain becomes closed. The result shows' .that 
the decay of the electrostatic wave energy is governed by the sum of * 
two.' damping coefficients: (1) a linear Landau damping coefficient 
which is independent of the amplitude, and (2) a turbulent damping 
coefficient dependent on the amplitude. Correspondingly the oreieuuriu 
coefficient consists of a linear part' arid a turbulent part. On the same- 
basis a kinetic equation is derived in the generalized form of the 

-x 

Fokker-Planck equation. 
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PLASMA VORTICES AND THEIR MOTION IN INHOMOGENEOUS MAGNETIC FIELDS 

Winston H. Bostick 
Stevens Institute of Technology- 

In the laboratory experiment of plasma flow over a magnetic dipole the 
breakup of the plasma front into plasma vortices has been observed by magnetic 
and electric field probes and by image-converter photography. (See figure l\ 
Recent image-converter photographs, now show quite clearly the formation of 
these columnar type vortices in pairs. One can make a plausible model for 
the formation of these pairs out of the flutes which grow in Rayleigh Taylor 
instability. (See figure 2), One ban also now understand on the same basis 
why the plasma proceeding across a magnetic field from a small two-wire button 
plasma gun is in the form of columnar plasma vortices two abreast. 

The photographed trajectories of plasma from a two— wire button plasma gun 
fired in inhomogeneous magnetic fields shows clearly that the plasma vortices 
"bounce off" of the region of increasing field. 

The force-free type of vortices with helical magnetic fields observed by 
Wells have been studied by one of our students Robert Small, at Grumman. • He 
confirms with density-flow probes, magnetic and electric field probes that the 
conical theta-pinch generates a train of such force-free vortices which proceed 
along a guide field and coalesce. with time. 

The plasma coaxial accelerator operated with a hexagonally -shaped center 
conductor generates two such force-free vortices on each flat side of the 
hexagon. We now understand qualitatively the process whereby these pairs of 
vortices, one corotating, the other counterrotating, are generated. (See figure 

3). 
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Addendum to "Report 

' > 

Pair Production of Plasma Vortices , W. H. Bostick et» al; 

The only vectors in figure 5 whose directions (i.e. the signs) have 
been experimentally measured or inferred thus far are: 

a. The general direction of B g in the coaxial accelerator, restating from 
a negative center conductor and positive outer conductor, 

b* The direction of B along each of the radial pl'asW filaments (see 
figure 3). 

c» The plasma density flow vector along the filaments is inferred from 

general plasma flow patterns in previous coaxial accelerators and 

also from the fact that figure 2 shows a piling up of the plasma at the' 
i 

locations where the radial filaments intersect the outer conductor (anode). 

In figure 5 the plasma density flow vectors circumferential to the filaments 
and magnetic fields circumferential to the filaments have "been only surmised with 
consistency being -observed" in righthandedness and lefthandedness for’ both magnetic 
and velocity vectors, A serious criticism of the directions chosen for these 
vectors in figure 5 is that the current along the filaments would be in the * 
opposite direction to the general current in the coaxial accelerator. An 
alternative and probably- preferrable choice is. that shown in figure 6 in this 
addendum. The true determination of the directions of these vectors must awain 
further measurements with coupling loops and density flow probes. These measure- 
ments are now in progress. 

ERBATA 

1, Page 2 line 12 should read B. in the coaxial accelerator) and show- that 
these these fields are of opposite — — — 

2. Page" 9 figure 5 caption - delete last line, 

Winston H* Bostick 
Stevens Institute of Technology 
Hoboken, hew Jersey 
5 April 19 66 




plasma density 
flow vector 
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EXPERIMENTAL UUVESTIGATION OF- THE FUIHDAMENTAL 
MODES OF A COLLISIONLESS PLASMA 

■Contract Humber NAS7-275 

General Atomic Division 

General Dynamics Corporation 

P. 0. Box 6083 San Diego ,* California 92112 

Principal Investigator : 'John H. Malmberg 

Associates: David L. Book 

Norris W. Carlson 
T. K. Fowler 
C . Dave Moore 
T. O’Neil 
Charles B. Wharton 
Norman Rostoker 

Objective 

To measure .the dispersion, damping, and heam induced growth of plasma 
waves near the electron cyclotron frequency, and to develop theoretical disper- 
sion r elat ions for plasma waves in non-uniform plasmas. This work is part of 
a project to investigate experimentally and theoretically the turbulence resul- 
ting from the growth of waves in a collisionless plasma, and to study the 
effect of the turbulence on properties of the plasma (e.g. , transport coef- 
ficients ) . 

Background 

The behavior of plasmas is, in many circumstances , dominated hy dynamic 
effects' which do not depend on the ' existence of large-angle collisions. In 
these collisionless plasmas the usual collisional dissipative mechanisms are 
absent and the dynamics of such plasmas are completely different from the 
dynamics of a fluid. Among the most basic properties of a collisionless plasma-, 
from a theoretical point of view, are those embodied in the dispersion relations 
which determine the character of the waves which can propagate in the plasma. 

If any of these waves are unstable, they grow to some limit determined by the 
non-linear theory of plasma dynamics, and the wave spectrum in the final state 
(i.e., the turbulent end-result of the growth) determines many of the plasma 
properties of practical interest. For’ example, the transport properties, the 
microwave scattering cross section, and the equilibrium density in a contain- 
ing field (e.g., the earth *s magnetic field) are likely to he dominated by the 
turbulence of the plasma. 

The dynamics of a collisionless plasma is described by the collisionless 
Boltzmann equation combined with Maxwell’s equations. These equations, com- 
bined with the definition of charge and current, constitute a complete set of 
equations and thus, if boundary conditions are specified, are in principle a 
complete description of the collective motions of the plasma. However, the 
equations are - non-linear and cannot he integrated directly. To obtain solutions 
the set of equations .is usually linearised by investigating the effect of a 
small disturbance on some initial state, i.e., using perturbation theory. By 
carrying the perturbation theory to second order ( "quasi -linear theory ) a 
description of the turbulence resulting from the wave growth is obtained, and 
from what amounts to a third order calculation, the ultimate dissipation of 



the turbulence is predicted. Experimental confirmation of the theory is al- 
most completely lacking, partly because it has only recently been possible to 
make laboratory plasmas "which match the assumptions of the theory. 

Progress and Results 

The principal theoretical effort in the program has been devoted to de- 
veloping the theoretical techniques necessary for 'the application of the linear 
theory of plasma waves to the geometries encountered experimentally, -where the 
size is finite, and the density a function of position. In particular, methods 
have been developed for predicting the dispersion ;and damping of electron 
cyclotron waves in. the geometry of the turbulence machine from first principles . 

The turbulence machine at General Atomic (developed under .other programs ) 
provides 1 a quiescent, collisionless plasma which matches the assumptions of the 
theory. Included in the installation are a variety of highly developed diagnos- 
tic instruments for wave propagation experiments. Many of the plasma properties 
are known in detail from other experiments. Using this, facility, the propaga- 
tion of waves near the. electron cyclotron frequency has been investigated. The 
dispersion and damping of the waves has been measured, and their growth due to 
an electron beam injected into the plasma has been observed. Some of the 
waves are wave-guide modes perturbed by the plasma. The wave of most Interest 
for these measurements is electron cyclotron wave. This is a backward Wave 
which propagates between the cyclotron frequency and" the upper hybrid frequency. 
The measurements have been made for a variety of .plasma parameters. 

Addition of a cusp magnetic field to the usually uniform magnetic field 
near the plasm source reduces the anomalous diffusion of the plasma by an 
order of magnitude. We hypothesize that this geometry decreases random electric 
fields in the plasma by reducing communication with the .source, and thus re- 
sults in reduced diffusion. 

Future Plans 

We intend to continue the investigation of the waves. We want to com- 
pare the observed dispersion relations with theory using independently measured 
plasma properties in the theory for comparison. The damping measurements 
need better precision, and much better growth data would be desirable. After 
the dispersion, damping, and linear growth of these waves is measured and 
compared to theory, the next step along the road to. understanding the ultimate 
turbulence and diffusion associated with them is to ■ investigate the non-linear 
limit of their growth. Further measurements will also be made to establish 
whether our hypothesis on how the cusp field reduces diffusion is correct. 
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STABILITY OF PLASMAS SUBJECT TO CONVECTIVE INSTABILITIES 


E. H. Holt and B. A. Huchital 
Rensselaer Polytechnic Institute ,, Troy, N.Y- 


A transverse magnetic field "which is "weak relative to the main longi- 
tudinal magnetic field applied to a positive column plasma has the effect of 
increasing the range of plasma parameters over "which the plasma remains stable 
against the well known Kadomtsev or convective instability. Kadomtsev and 
Redo spas ov- have shown that a " perturbation of the form 

n(r) = J-^Sr) exp (im$ + Iks + lot) 

will grow provided the criterion 

Ak 4 +' Ck 2 + G + C < mHku E XT (l) 

o tn s z © 

is satisfied., where the coefficients are functions- of the longitudinal magnetic 
field defined in reference 1. ¥e have modified equation (l) to include a 
transverse field, by writing it in the form 2 

AkVp 2 + Ck 2 /P + G + G ^ mHku E /PT (2) 

' 'om e s' e 

where P - 1 + (p B ). 2 Clearly equation (2) reduces to equation (l) in the 
limit B^O. 

The two criteria are compared in Figure 1, in terms of the parameter 
9 = tan B /B . It can be seen that the plasma becomes more stable as B 

increases. Experimental confirmation of this result has been obtained. A 
further consequence of finite 2L is shown to be the surpression of m = 1 
in favor of m - 2 instabilities . 3 

Results of experiments aimed at stabilizing the plasma by means of 
quadrupole magnetic fields are illustrated in Figure 2. Stabilization of solid 
state plasmas in XnSb has been successfully achieved in this way by Ancker- 
Johnson. • In our ease, modest currents through the stabilising windings are 
shown to quench the instability oscillations quite effectively. 

The convective type of instability is likely to occur in a number of 
different plasma types covering a wide range of plasma parameters. Guest 
and Simon5 have suggested that an entirely similar mechanism operates in the 
low pressure arc due -to streaming of particles along the magnetic field. 

Holi and Simon? have proposed that a variation of the Kadomtsev instability 
is present in a reflex discharge. Finally, Cherrington^ has shown that 
longitudinal gradients of density and temperature pan trigger the Kadomtsev 
instability, thus suggesting an explanation for the anomalies observed in an 
rf discharge hy Powers.? A systematic approach, based upon the techniques 
we have developed in studying the positive column, is expected to yield 
important information concerning the instability mechanism in this wide variety 
of plasmas. 
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Figure 1 . The Stability Criterion for a Positive Column Showing the Effect 

of a Transverse Component of the Magnetic Field. 

(tan 9 = B. /b, ) 

transverse 7 longitudinal 7 



Figure 2. 
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MODES OF THE KADOMTSEV INSTABILITY 


D. A. Huchital and E. H. Holt 
Rensselaer Polytechnic Institute, Troy, N.Y. 


We have carried out a systematic experimental study on the modes of 
the Kadomtsev instability of a collision dominated magnetoplasma which has 
resulted in new information and understanding concerning the presence and 
dominance of the separate modes, m = 0, m = 1 and m = 2. 


Anomalous effects were reported by Lehnert to appear when a longitudinal 
magnetic field, B^, impressed on a dc positive column exceeded a critical 
value, B c * The onset of these effects was explained by the stability con- 
siderations of Kadomtsev and Nedospasov 2 who showed that a helical perturbation 
of the form (f3r) exp(im0 + ikz + icnt) will grow exponentially when B > B 
and m 1. A L c 


We have studied the plasma with Langmuir probes to determine the rela- 
tive importance of the various possible modes. 3 We have used sets of four 
probes inserted 90° apart to measure azimuthal phase relations. A typical 
result is shown in Figure 1. The 90° phase difference between the two signals 
verifies the presence of the m = 1 mode. 

Our investigations have included the interesting observation that an 
m = 0 class of oscillation is commonly present in the plasma. This new result 
cannot be explained on the basis of the Kadomtsev theory. Data will be 
presented to show that the m = 0 oscillation vanishes when the electric 
field in the column is precisely aligned with the main magnetic field. 
Transverse components of magnetic field as small as 0*5 gauss are sufficient 
to generate m = 0 oscillations of appreciable amplitude. Experimental and 
theoretical results will be presented to show that the m = 0 oscillation is 
an important non-linear effect resulting from the interaction of the helical 
m = 1 mode with the transverse magnetic field. 

When the transverse magnetic field is increased to the order of several 
tens of gauss, the m = 2 mode is observed to become the dominant oscillation. 
However, the experiment is incapable of distinguishing between a genuine 
m = 2 instability and a non-linear effect similar to that associated with 
m = 0. We have developed a new experimental technique^ that considers the 
response of the plasma to externally impressed perturbations while the magneto- 
plasma is operating in the sub-critical (B^ < B c ) regime. 

Figure 2 shows the signal detected at probe set B as a function of 
the frequency of the source connected to probe set A. It is shown that as 
the critical magnetic field of 2100 gauss is approached a distinct resonance 
representative at the m = 1 instability is excited. Experimental results 
obtained by this technique will be presented to establish the independence 
of the m = 2 instability and to confirm the non-linear nature of the m = 0 
mode. 
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figure 1 . Probe Signals Illustrating the 
Helical Mode of the Instability. Top 
trace: V ' - V . Bottom trace: 

Time scate: 20 ps/cm. • Figure 2. Signal -Received at Probe 

Set B as a Function of the Frequency 
of the Source Connected Across Probe 
Set A. At the critical magnetic 
field strength (2100 gauss), the plasma 
exhibited seif sustaining oscillations 
at 16. 3 kilocycles. 




PLASMA RADIATION AND THE DETECTION OP NON-MAXfflYLLIAN 
ELECTRON VELOCITY DISTRIBUTIONS* 

J.H. Noon, P.R. Blazuk, E-H. Holt 
Rensselaer Polytechnic Institute, Troy, N.Y. 


A gated, microwave radiometer and a specially constructed magnetoplasma 
■waveguide cell^ have "been used in this laboratory to monitor the radiation 
temperature of electrons in an afterglow plasma and to examine departures 
from a Maxwe lli an electron velocity distribution. The cell is also suitable 
for a study of anomalous diffusion in a pulsed magnetoplasma, or for measur- 
ing noise radiation near the critical magnetic field for plasma instability, 
although these measurements have not as yet been carried out. 

Microwave radiometers have "been used by a number of workers 1 to detect 
the thermal noise radiation from electrons in a collision-dominated nitrogen 
plasma and thus to determine the electron temperature, assuming a Maxwellian _ 
electron velocity distribution exists in the plasma. This assumption is not 
necessarily valid always as shown by our results which are reported below. If 
a - longitudinal magnetic field is applied, measurement of the noise power from 
the plasma, as a function of magnetic field, for a range of field strength 
such that the electron cyclotron frequency is close to the microwave observa- 
tional frequency, may be analyzed^ to determine departures from a Maxwellian 
distribution. A sample of our experimental results in a nitrogen discharge 
showing noise radiation temperature as a function of applied magnetic field, 
are shown in Figure 1 , and distribution function parameters giving best fit 
to the data are shown in Figure 2 - These results show clearly that for a 
short time after the cessation of the pulsed DC discharge the nitrogen after- 
glow is relatively rich in high energy electrons, then becomes deficient in 
high energy electrons, then relaxes back to a Maxwellian distribution. This 
also manifests itself as an apparent maximum in the radiation temperature, at 
post-discharge times of order 25 microseconds, and by a relatively slow rate 
at which the electrons relax back to room temperature. The magnitude of both 
these effects depends somewhat on the energy input into the discharge as shown 
in Figure 3 * Time resolved studies of the active discharge itself have also 
been carried out and it appears that a time of the order of 10 microseconds 
after application of the DC voltage to the electrodes is required before a 
Maxwellian electron velocity distribution exists in the nitrogen discharge. 

These, results can be used both to provide information about energy loss 
processes in the plasma and also to indicate the regime in which the magneto- 
plasma becomes unstable above a critical magnetic field. 

Further measurements planned with the magnetoplasma cell will be outlined. 


*Work supported by NASA under grant no. NsG 48 . 
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Figure. 1. Radiation noise temperature as a function of applied magnetic 

field B. When A(B) = 0, B is such that the electron cyclotron 
frequency cn^ is equal to the microwave observational frequency co. 

H D p p 

A(b) = 10 ((cn - o^)/P Q ) (cycles/sec torr) 

Figure 2. ‘Variation of distribution function parameter^ with time in the 
afterglow. Assuming f (v) oC exp(-bv^J, JL - 2 corresponds to a 
Maxwellian and Xj = 4 to a Druyvestyn distribution. 

Figure 3* Radiation noise temperature in the afterglow period for different 
energy input conditions In the active discharge. 

X= O. 5 A, • = 5A, O = 7 A, and all are 20 microsecond pulses. 
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PLASMA MEASUREMENTS PROM 3 TO- 120 KILOBARS OF PRESSURE 

James V. Robinson 

The University of Michigan, Ann Arbor, Michigan 
Grant No. NsG-4l5 


1 

The purposes of the research have been to make basic measurements 
of the properties of very dense plasmas and to check the accuracy of 
theories "which might be applicable to a description of these plasmas. 
Measurements of temperature, pressure, resistivity, and energy density 
were obtained in the range of pressures from 3 to 120 kilobars. The 
Debye shielding theory often used for theoretical studies of plasma was 
found to be very inadequate at these high 'pressures. 

Although high density plasmas occur in various experiments involv- 
ing exploding wires and hypervelocity guns, for example, little has 
been known about the properties of such plasmas. E. A. .Martin of the 
University of Michigan introduced techniques which have been extended 
to form the basis of the present research into the properties of plasmas 
at high pressure. 

The experimental work was conducted in two phases. First, with 
a refinement of Martin's methods, plasmas were studied in the range of 
pressures from 3 to 15 kilobars. These plasmas were formed by discharg- 
ing a capacitor bank between electrodes submerged in water. From mea- 
surements of radiation intensity, measurements of the electrical dis- 
charge characteristics, and photographs, the desired data was computed. 
The pressure was a function of the rate of current rise in the column 
of pi asms, formed by the discharge. In the second phase, the discharge 
was formed in water which had been compressed by a chemical explosive 
to 100 kilobars of pressure . The data was recorded similarly and varia- 
tions of pressure from 100 to 120 kilobars also depended upon the rate 
of current rise. 

The diameter of the discharge column and the current flowing 
through the col umn increase nearly linearly with time . As the column 
grows, energy Ib supplied from the discharge circuit so as to maintain 
a nearly constant energy density In the plasma. Likewise the other in- 
tensive properties of the plasma are nearly constant with time. Energy 
lost from the plasma through expansion against its surroundings ^ and , 
through radiation is a small fraction of the total energy supplied to 
the plasma. The best data are taken late in the interval of observa- 
tion when the variables are in the most accurate ranges of the instru- 
mentation, that is in the interval from 0,8 to 1.0 psec. 

The magnetic pinch effect is important In the calculation of pres- 
sure. For a unif orm radial current distribution the pinch pressure 
varies from zero at the boundary of the discharge to a. peak on the cylin- 
drical axis. The volume average of the pinch pressure which is half of 
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the peat: is added to pressure produced by the inertial restraint of the 
water surrounding the column to obtain total pressure. Since the resis- 
tivity of the plasma has been found to decrease with increasing pressure, 
the radial pressure gradient tends to concentrate current at the center, 
causing an even greater pinch than has been estimated. 

Without exception, the temperatures of the plasmas formed without 
explosives were about 35 000°K while the temperatures with explosives 
were about 10 QQ0°K. The variations with pressure of the energy density 
u and the resistivity p are plotted in Jig. 1 and Fig. 2 : 


A system of equations based upon the Debye shielding theory was 
used to predict particle density and energy density from temperature 
and pressure for the case without explosive. The predi'cltiong of energy 
density were found to be about one -fifth of the measured values, while 
the measurements are considered accurate to within 20 percent. The 
failure of the Debye theory was not surprising because the Debye radius 
was smaller than the average interparticle spacing. In place of the 
Debye theory, a semiquant itative theory based upon the distortion of 
electron quantum levels was found to provide a much better explanation 
of the energy density measurements. For the case -with explosive, the 
Debye theory predicted an unreasonably low level of ionization because 
of the low temperature, and the analysis was not carried further. The 
concept of excluded volume was introduced into calculations for the 
higher pressure . From the theoretical considerations, the density of 
atoms for a 100 kilobar plasma was est ima ted to be 10^5 e m"5 and for a 
7 kilobar plasma to be 10 cm - - 3 . 


Publications 


Robinson, James W. 1 Measurement and Interpretation of -Plasma Properties 
to 100 Kilobars of Pressure, To be published as a NASA Contractor 
Report, CR~Ui|.6. 

Robinson, James W. : Plasma Measurements at 100 Kilobars- of Pressure. 

Appl. Phys. letters, vol. 8, no. 8> I3 April 1966. 
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N INVESTIGATIONS ON THE MECHANISM OF THE MATERIAL RELEASE 
IN THE HIGH VACUUM BREAKDOWN 

by 

Richard T, Schneider 

Department of Nuclear Engineering Sciences 
University of Florida 


ABSTRACT 

The program to be reported started on May 1, 1966; therefore, 
it is necessarily still in the planning stage. However, preliminary 
experiments for orientation purposes have been made. The program is 
aimed to make a contribution to the understanding of the mechanism 
of the high voltage vacuum breakdown. Special emphasis is placed on 
the mechanism of material release from the electrodes. It is desir- 
able to find means to predetermine the material release, and perhaps ' 
initiate it by magnetic ignition. The information obtained should 
he useful for electric propulsion using metal plasmas. 

The measurements will be made in a vacuum as low as 10”-^ 

Torr, The electrode configuration will be parallel plates for 
electrode conditioning and repeatability tests. Other configura- 
tions are planned to favor or suppress one or more of the many pro- 
cesses which take part in the vacuum breakdown. These processes 
are (1) electron frees positive ion at the anode, (2) photon frees 
electron on the cathode, (3) electron frees photon on the anode, 

(4) positive ion frees electron at the cathode, (5) negative ion 
frees positive ion at the anode. Using appropriate electrode con- 
figurations , with auxiliary electrostatic or magnetic fields , it 
should he possible to suppress or enhance some of these processes 
in order to study their influence on the mechanism of the vacuum 
breakdown. 

Spectroscopic studies of the inter- electrode plasma are 
planned. The object is to determine the temperature, density and 
ionization degree of the material released from the electrode. In 
order to do this, methods of short-time spectroscopy have to he 
applied. The discharge time is in the order of lysec. or shorter. 
Therefore, it is worthwhile to investigate whether the inter- electrode 
plasma has a defined temperature. Our point of interest will be to 
find out if the material released from the electrodes is released in 
the form of ions or if it is released as atom vapor by thermal vapora- 
tion, or if it is released in larger units (clumps) by mechanical 
stress or impact. 

Time integrated spectroscopy can be used to determine how 
much material from which electrode is released. This will be aided 
by making anode and cathode of different materials. 
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Some results of the preliminary experiments will be shown 
at the meeting. 



53 


SYNTHESIS AMD CHARACTERIZATI ON 
OF MAGNETIC FLUIDS 

R.E. Rosensweig 
Avco Corporation 

Research and. Technology Laboratories 
Wilmington, Massachusetts 


In the present research program ve endeavor to synthesise fluids 
having a strong magnetic response, and to characterize their properties 
and discover the laws and relationships which govern their response to 
an applied magnetic field. 

Background Information 

All known materials that are strongly, magnetic like "iron — the 
ferromagnetic metals, alloys and compounds — are solids. While theore- 
ticians do not rule out the possibility of true ferromagnetic liquids, 
none is known at present. It is, however, possible to make a stable 
suspension of a finely divided ferromagnetic solid in a carrier liquid, 
to produce a fluid that reacts strongly when a magnetic field is applied 
to it, and we are engaged in studies of such fluids. 

The pressure distribution in these fluids can be independently 
controlled by an applied magnetic field. The fluid response then takes 
several basic forms', opening up a wide range of possibilities in viscous 
dampers, inertial sensors, heat engines, fluid computers, fluid 
modelling studies, medical technology, and other technological arcaD.. 

Investigation of the basic fluid mechanical processes which 
govern the response of a magnetizable fluid to a source of magnetic 
field leads to an augmented Bernoulli principle in which a magnetic 
pressure-like term is additive to the usual pressure, speed, and 
gravitational terms . 

H 

p + •§■ q 2 + ^ gh - J* IdH = constant 

0 

The various terms in this equation along with 'consistent c.g.s. units 
are, p = pressure (dynes/cm^) , density, assumed .constant (gm./c.c.), 
q = speed (cm. /sec.), h = elevation (cm.), I = magnetic moment per unit 
volume (pole-cm.) = M/^TT where M is expressed in gauss (M = B~H), and 
H = applied field (oersteds). The initial three terms taken in various 
combinations describe fluid mechanically a number of familiar techno- 
logical processes as categorized in an accompanying figure. The 
magnetic . term then combines to yield additional classes of fluid 
dynamic and static phenomena laden with possibilities for technological 
application. These furnish a portion of the motivation for these 
studies . 
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Program Activities 

Mach of the work under the program has concentrated on the 
preparation and study of magnetic fluids composed of ferrite solid 
particles dispersed in an organic liquid. The carrier fluids are 
chosen on the "basis of viscosity, chemical nature, and "boiling point. 
Dispersants are selected on the theory that a suitable molecular 
structure contains a chainlike nonpolar tail and a polar end group. 

Both magnetite and a lower Curie temperature, manganese zinc ferrite 
are used as the solid, magnetic constituent. Dispersion into colloidal 
form is accomplished "by wet hall milling of the -solid material in the 
‘presence of the carrier liquid and the dispersant. 

Screening tests are conducted to determine systems that lead to 
the successful production of colloid. Then study is made of the effects 
of initial charge and the grinding time on the formation of colloid and 
its physical nature. The successful fluids are characterized by their 
properties of magnetism, viscosity, temperature dependence of viscosity, 
density, stability, and particle size and shape. The fluids are 
modified by vacuum distillation or addition of solvents and study is 
made of the modified fluids so obtained. 

Correlations are studied of the inter-relationship between the 
various measured properties and known parameters such as dispersant 
chain length and attempt is made to feed back this information into the 
synthesis of additional fluids having accented properties. 

Presentation of Results 

Stable magnetic fluids are produced in which the typical parti- 
cle size is 100 angstroms. The number of colloidal particles per cubic 
centimeter is enormous, and characteristically ranges from 10 1 ? to 10^- . 
Oleic and linoleie acids have proven successful as dispersing agents and 
so has another class of dispersants consisting of succinic acid deriva- 
tives . Fluids prepared in this manner are strongly attracted to a 
magnet and volumes of the fluid may be suspended against the force of 
gravity. The magnetic particles experience a constant attraction in 
the direction of an applied field gradient but at the same time are 
free to move through the fluid, hence the persistence of a stable sus- 
pension requires explanation. It is believed that at equilibrium the 
terminal motion of the particles equals the diffusional transport of 
particles due to a (small) concentration gradient. Then the stable 
diameter d is less than the value computed from 

d3 = 2k 

where k = Boltzmann's constant = 1.38 x 10 "^ erg/part iele °K, T = 
absolute temperature in Kelvin unit, M = ferric induction in gauss, 

H = applied field in oersted, and s is length — d is then given, in 
centimeter’s . 

Suspended in the fluid each particle is analogous to a molecule 
of a paramagnetic gas. The Langevin equation of statistical mechanics 
then relates the average magnetic moment per particle averaged over 
all thermallv induced orientations to 1) the true magnetic moment of a 





particle, H the applied field, and T the temperature. Unlike a para- 
magnetic gas molecule, for a particle of the magnetic fluid \) = >)(t) « 
M extended form of the equation has been obtained to deal -with the 
actual particle size distribution function. 

Hie experimentally determined dependence of riscosity upon the 
loading is correlated -with theoretical predictions . The earliest 
theory was derived hy Einstein who solved the flow field of pure 
strain as perturbed hy the presence of a sphere . His result relates 
mixture .viscosity 3^ to solvent viscosity and solids fraction <f> : 

|=l +2 . 5 ^. 

However, this relationship is valid only for dilute colloidal mixtures 
a better expression, correct for high solids loadings, has the form: 


3»- i 

% 1-2.5 j> + l-55^ 2 


This relationship predicts solidification corresponding to spheres in 
the hexagonal close-packed arrangement. In using this relationship, 
the thickness S of the monolayer coating on each particle must he con 
sidered in the computation of . „ 

Another important consideration is the prevention of particle 
agglomeration or flocculation. The important mechanisms are the mag- 
netic forces and van der Waal’s attractions (which tend to flocculate.) 
and repulsion forces arising from deformation of an adsorbed layer of 
dispersing agent. The energies of two adjacent particles are found to 
be: ■ 


Magnetic energy: 


Van der Waal's: 
(London's' model) 
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Where 1'T is the number .of adsorbed molecules per unit area, S is the 
length of the hydrocarbon chain, and h is the relative separation of 
particle surfaces, h = H/r where H = surface to surface’ separation. 

For stability against flocculation the three energies must 
produce a superimposed energy distribution having a harrier or energy 
hump at least several times greater than KD, the thermal energy of a 
particle . 
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At the beginning of this program, magnetic fluid with, viscosity 
less than ten centipoise was limited to a ferric induction at satu- 
ration of little more than 100 gauss, The relationship between vis- 
cosity and ferric induction is presented lit the accompanying figure 
which also illustrates the advance in the state of the art. Fluid is 
now produced whose ferric induction exceeds 1000 gauss . The remaining 
figure presents data that permits correlation of the initial perme- 
ability of the fluids with the ferric induction at saturation, Mg, 
according to 


A- ^ 1 + 


Ms N 

mo } 


Mg > 100 gauss 


This tells that the present fluids yield permeabilities in. the range 
1 </it< 3 which, although it is extraordinarily large compared ,to the 
permeabilities of normal liquids, is still rather small compared to the 
permeabilities of solid ferromagnetic materials . The obtaining of 
stable magnetic fluids having a much greater permeability is one of a 
number of challenging problems in this work where the solution in 
probability depends upon an understanding of surface physics and 
chemistry, magnetostatics, and the behavior of matter in general. 


Background References 

Rosensweig, R.E., Nestor, J.W., and Timmins, R.S., "Ferrohydrodynamic 
Fluids for Direct Conversion of Heat Energy’* A.I.Ch.E.-I . Chem. E. 
Symposium Series No, 5? London, p. 104 ( 1965 ). 

Neuringer, J.L. and Rosensweig, R.E., "Ferrohydrodynamics”' The Physics 
of Fluids 7 No. 12, p. 1927* December (1964), 

Rosensweig, R.E., "The Fascinating Magnetic Fluids" New Scientist p.l46 
January 20, 19 66 . 

Publications Generated During This Research 

Rosensweig, R.B., "Buoyancy and Stable levitation of a Magnetic Body 
Immersed -in a Magnetizable Fluid." Accepted for publication in Nature , 
A Weekly Journal of Science. 

Neuringer, J.L., "Magnetic Effects on the Diffusion of Ferromagnetic' 
Colloidal Particles Dispersed in a Current-Carrying Fluid Medium. ” 

To appear in The Physical Review Part I in May ( 1966) - 
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PROPAGATION- AND DISPERSION OF HYDRO MAGNETIC AND 
ION CYCLOTRON WAVES IN PLASMAS IMMERSED IN MAGNETIC FIELDS* 
Professor Arwin A. Dougal, The University of Texas /Austin 78712 


The long range challenges for exploration and exploitation of space, requires 
that ever increasing, technological and scientific attention be given towards the 
realization of extremely light, almost infinitely long-lived, and efficient power 
sources. Within, very recent years the technological community has recognized 
the vast potential power production realizable from the fusion products of reacting 
light nuclei as the ultimate source. Almost overwhelming power output per unit 
fuel mass is afforded by the reacting light nuclei. Two million kilowatt hours- of 
energy are forseen to be extractable through fusion of as little as 0,06 pounds of 
deuterium nuclei. An immense technological problem to be overcome is the require- 
ment for the production and heating of deuterium ions in the gaseous plasma state 
to several hundred million degrees temperature to which this research work and its 
continuation are addressed. Experimental and analytical findings from our recent 
research work on ion cyclotron resonance heating show definite promise fowards- 
future realization of increasingly super-high ion temperatures , These findings are 
supported by relatively recent reports from those few research groups throughout 
the world who have'had the foresight to initiate studies on ion cyclotron resonance 
heating and to- realize the major experimental facility for effective experimental 'work. 

Most earlier and existing attempts to elevate plasma temperatures into the 
desired range have resulted in high electron temperatures but relatively low ion 
temperatures which are inadequate - Both theory and experiment have recently con- 
firmed' that effective heating, propagation, and dispersion of- ion cyclotron waves 
in plasmas immersed in magnetic fields are possible. Deposition of power into a 
plasma through radio frequency fields at an angular frequency CD requires con- 
sideration of characteristic frequencies associated with the plasma and its magnetic 
field configuration. The characteristic frequency ~C2-jL of an ion in a magnetic 
field is given by At » , where q i and nr are the ionic charge and 

mass respectively. Theory and experiments show that effective heating of plasma 
ions occurs through ion cyclotron wave interactions , 

Experimental research on ICRH has been performed by research groups at NASA- 
Lewis , Piincetbn University, Lawrence Radiation Laboratory at the University of Cal- 
ifornia at Berkeley, Physicotechnical Institute at Khar'kov USSR, The University of 
Texas, the Massachusetts Institute of Technology, and the Institute of Plasma Physics 
at Nagoya University, Japan. 

The experimental arrangement at The University of Texas is shown in Fig. 1, 

Both 'hydrogen and/or deuterium gas can be leaked into the plasma tube; The vacuum' 


*Sponsored by the National Aeronautics and Space Administration Grant NsG-353. 
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base pressure is better than 10” torr. VacSorb and Vaclon pumps are used in 
order to reduce impurities in the system. The overall length of the magnetic field 
coil system is presently 150 cm. Rf power is coupled into the plasma through a 
Faraday shielded "Stix" coil located in a uniform field region of maximum 10.6 kG. 

A "magnetic beach" region on one side of the "Stix" coil is used for wave propa- 
gation studies. Maximum incident rf power is presently over 30 kW. The incident 
and reflected rf power is measured with adjustable, calibrated directional couplers. 
The rf driver is crystal controlled at 5,8 MHz. A 27 MHz, 1 kW, preionizer is 
used in all experiments. A fairly extensive study of the design of the "Stix" coil 
and associated electrostatic shields was undertaken. The present coil and shield 
combination gives an E$ field which is almost sinusoidally varying in the axial 
direction and an field which is substantially reduced. The number of coil turns 
and their relative spacing was determined using a digital computer. The Faraday 
shield consists of 9 strips of copper, .020" thick by 0,5" wide along the discharge 
tube under the Stix coil. 

A new and previously unreported wave phenomenon was observed in this ex- 
perimental investigation of ion cyclotron resonance heating of hydrogen and deute- 
rium plasmas. For certain gas pressures the ion 'cyclotron waves propagate for 
frequency (f) greater than the ion cyclotron frequency (fj ) , and attenuate at har- 
monics of the ion cyclotron frequency with a resulting increase in the transverse 
temperature-density product (nKTj ). In. this experimental investigation the fre- 
quency (f) of the driver was kept constant and the other parameters , such as gas 
pressure and magnetic field strength, were varied. Decreasing the magnetic field 
strength is effectively the same as increasing the frequency of the driver (f) , and 
vice versa. The wave magnetic field components were measured with movable, 
shielded, magnetic probes located near the end of the magnetic beach. 

For low hydrogen pressures (p< 2 millitorr) an ion cyclotron wave was observed 
to propagate for B ->B C , where B c is the ion cyclotron resonance magnetic field for 
the driver frequency (f) . The ion cyclotron wave was observed to attenuate quite 
sharply near B = B c as expected. As the gas pressure was increased, a wave was 
observed to propagate for B < B c , which was unexpected. Closer investigation re-, 
vealed that the wave attenuates near B c /iT> where N is an integer. Thus, the wave 
propagates for f (the ion cyclotron frequency) and attenuates at harmonics of f£. 
The value of N was found to depend upon the gas pressure such that N- increased 
for higher pressures. The same general effects were observed in deuterium plasmas. 
In deuterium, the effects are more pronounced since the separations between the 
various B c /N are larger due to the larger deuterium ion mass. Fig. 2 shows these 
effects as measured in a deuterium plasma. The region of abrupt wave attenuation 
moves to lower magnetic fields as the gas pressure increases. The attenuation 
regions are very close to the harmonic resonance fields . 

The question arises whether the observed wave field attenuation of harmonics 
of % is. due to a wave cut off or to actual wave damping near the critical fields , Dia- 
magnetic probe measurements indicate that the waves are damped and that energy is 
transferred from the waves to the plasma particles . Additional measurements on 
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wave propagation and wave attenuation are in progress. An analytical explan- 
ation for these observed effects is being developed. 
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A-C POWER GENERATION THROUGH WAVE -TYPE INTERACTIONS 

by 

G. L. Wilson 

Department qf Electrical Engineering 
Massachusetts Institute of Technology 


The purpose of this project is to investigate schemes which 
directly generate A-C electric power through magnetohydrodynamic 
interactions in plasmas. The method under current study is one 
in which standing magnetoacoustic waves -in a flowing plasma inter- 
act with an artificial or lumped transmission line, ■sdiich is 
magnetically coupled to the gas. Magnetoacoustic waves are distur- 
bances which propagate in a conducting gas at velocities greater 
than the sound speed in the medium. This greater propagation velo- - 
city is due to currents which are induced through the compression 
of magnetic field lines within the conductor. These currents 
interact with the applied magnetic field to produce forces which 
aid compressional forces. 

The proposed generator is shown in Figure 1, in which the 
coils of an artificial transmission line are coupled to the mag- 
netic field B associated with the waves in a convecting plasma. 

For a velocity V greater than the magnetoacoustic velocity, there 
will be. two forward waves in the gas. By adjusting the trans- 
mission line parameters such that the phase velocity on the trans- 
mission line is equal to the velocity of the slow wave in the gas, 
an interaction between the two systems occurs, which results. m 
unstable bulk oscillations in the plasma and unstable ^ voltages and 
currents on the transmission line. By coupling to this system 
in a variety of ways, net electric power can be coupled out of the 
system at the expense of the thermal and kinetic energy of the, 

' plasma. • 

Experiments are in progress in which a helium plasma is heated 
and accelerated in a circular channel through electric heating an 
T x B forces. The gas is accelerated to velocities on the order 
of 20,000 meters/second and is heated to temperatures in excess ot 
18 000°Kelvin. The experiment lasts for one millisecond at input 
powers on the order of 3 megawatts . A two-dimensional Hartmann 
flow theory for the accelerator has been developed and compares 
very favorably with measured values of the accelerator parameters. 
The results of these measurements are presented in Figure 2. 
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Measurements, of gas temperature and velocity profiles are in pro- 
gress to further check the Hartmann flow model.. 

The accelerator is used to provide proper gas conditions for 
study of magnetoacoustic wave oscillations. Magnetoacoustic waves 
have been successfully excited and detected in this medium and the 
acoustic -magnetoacoustic transition region as a function of elec- 
tric conductivity has been observed. The measured frequencies of 
the standing whve oscillations as a function of input power to the 
accelerator is shown in figure 2, and compared with theoretical 
predictions. This experiment represents the first instance in which 
magnetoacoustic waves in a flowing medium and- their dependence on 
electrical conductivity have been observed. 

Preparations are presently in progress for the final phase of 
the study; that is to experimentally observe the interaction of an 
artificial transmission line with the waves and to extract electric 
power from the system. 
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ANODE OSCILLATIONS .IN CESIUM PLASMAS 
Harry S. Robertson, Department of Physics, 
University of Miami, Coral Gables, Florida 


Self-excited oscillations in otherwise quiescent cesium 
arcs have been found to originate at the anode surface. ' 

These oscillations appear to be similar to those observed 
in cesium -diodes and in Q-machines with ion-rich sheaths. 

The suggestion that the oscillations in cesium diodes might 
provide a means for direct conversion of thermal to a.c. 
power on the one hand, and the experimental evidence that 
some of the instabilities seen in Q-machines originate in 
the anode sheath have motivated the study of these oscilla- 
tions in an experimental facility that permits a variety 
of observations. 

The experimental facility consists of a pyrex discharge 
tube with {usually) a thermionic cathode, usually about 30cm 
from, one or more of a variety of anode structures. A cesium 
reservoir in its own oven determines the vapor pressure, and 
the discharge tube is operated in a transparent pyrex fur- 
nace at a temperature usually 50°C or more above that’ of the 
reservoir. The neutral cesium pressure may be varied from 
about I0“ 3 Torr to 1 Torr; the current from a few m-illi-^ 
amperes to 1* ampere; the electron temperature is a few tenths 
of 1 eV, and the fractional ionization in the plasma is 
usually less than 1%, Anodes may be plane, cylindrical, or 
spherical, and multiple anodes recessed in cylindrical pyrex 
tubes are ^ sometimes used. <Both rubidium and potassium have 
been studied., but the most extensive work has been in cesium. 

Oscillation amplitudes may be as high as several volts, 
and periods vary from- a few- microseconds to a few milli- 
seconds. Photoelectric observations of the emitted light 
Shows that ^ the oscillations are localized to the’ region with-' 
a ew millimeters of the anode surface. 

, Theoretical and experimental analyses -lead to the follow- 
ing qualitative- understanding of the phenomenon. Neutral 
atoms incident on the anode surface are ionized because the 
work function of the anode is greater than- the ionization 
potential of the. atom, but they are lightlyhbound to the 
anode surface by image forces. A rising electric field at 
t e anode results in a burst of emitted ions from the surface, 
accompanied by release of electrons from the sheath edge to 
maintain neutrality. The resulting enhanced conductivity of 
the anode sheath leads to an abrupt fall in tube volts and 



local electric field, while the locally increased ion oh 
electron densities lead to a sharp increase in the emission 
of both recombination radiation and excitation radiation. 

Then, as the burst of ions and electrons dissipates, prin- 
cipally by diffusion and conduction, but in part by recom- 
bination, the electric field begins to rise. During the low- 
field part of the cycle, more neutrals have been absorbed on 
the anode surface, and when the field increases sufficiently', 
a new cycle begins. 

The oscillations will exist, then, whenever the steady- 
state electric field is within a range such that it is great 
enough to cause some field emission of anode surface ions 
and not so great that ’it steadily removes them as fast as 
they are formed. Frequency is determined principally by the 
diffusive and convective loss rate. 

The greatest difficulty in comparing theory and experiment 
is that ion emission usually occurs at one or more spots on 
the surface where the , local electric field must be a "bit 
larger than anywhere else, so brightly glowing spots are seen 
on the anode, and the firing of one spot triggers that of 
others. 

Attempts are in progress to control the formation of spots 
and otherwise alter the behavior of these oscillations by 
choice of anode geometry, material, finish, temperature, and 
geometry of the^ surrounding region. 
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Transverse Stream Instabilities in. Plasmas 

Dr. W. Bennett 
North Carolina State Univ. 

The objective in this project is to produce a steadily running 
undriveh pinch and to examine the instabilities which may exist in such 
streams. 

The project is supported by the Electrophysics Branch of the 
National Aeronautics and Space Administration and involves the production 
of high current high voltage electric pinches in small diameter beams and 
to observe the transition of these beams into a. simultaneous electric and 
magnetic pinch. Studies are being made of: (l) the physical processes 

in these two kinds of pinch; (2) the onset of beam instabilities as the 
magnetic pinch becomes dominant; (3) the interaction of these beams with 
their associated plasmas. 

In earlier work, it was found that undriven -electron beams do 
not readily accumulate ions from low density residual gas and retain them 
in sufficient quantity to neutralize the space change and permit the 
magnetic pinch to set in. For this reason, it was decided to undertake to 
increase the intensity of electrically pinched beams sufficiently to provide 
the essential conditions for the magnetic pinch to become predominant. 

In order to undertake this kind of investigation properly, it has 
been necessary to develop a directional plasma probe with which velocity 
distributions can be measured. The basic principle used in this probe is 
similar 'to that in a probe which has been successfully flown on a rocket 
in some earlier space research. The probe in its present form consists of 
two knitted 0.0005 in. tungsten grids and a 15° conical collector. The 
outer grid is attached to a metal sleeve which surrounds the rest of the 
probe. The leads to the interior parts of the probe are glass covered. 

The outer grid and sleeve core are held at the potential at which the probe 
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least disturbs the surrounding plasma. For measurements of the electron 
velocity distribution in the direction of the probe axis, the inner grid 
is held at positive potential sufficient to turn back all positive ions 
which enter through the outer grid. The electron current to the conical 
collector is measured with a full feed hack amplifier based at potentials 
which are swept through a range of negative values to cover the range of 
electron velocities in that direction. Corrections must be made for sec- . 
ondary emissions from the collector and for currents due to ionization 
inside the probe. For measurements of the positive ion distributions, the 
inner grid is set at a negative potential sufficient to turn back all electrons 
which pass the outer grid, and the base potential of the amplifier is swept 
through a range of positive values to cover the range of ion velocities in 
the direction of the probe axis. 

A tube has been built for examining electrically pinched streams 
by measuring velocity distributions of electrons and ions with directional 

prohes, magnetic stream deflections, stream attenuation, and other character- 

\ 

istics as a function of vapor density, beam current, beam voltage and tube 
diameter. 

Methods have been developed in the gun design for protecting 
coated cathodes from deterioration due to back bombardment by positive ions. 

In essence, the procedure is to use a "grid" aperature as near the cathode 
as possible which is held at a more positive potential than any other poten- 
tial further down stream. In order to go to the high voltage streams which 
will be required in the later phases of this project, where much more 
positive potentials will be necessary beyond the end of the gun, it will be 
necessary to design the gun to strongly defocus the back- streaming positive 
ions before they cair^reach the coated cathode, or else to use bare cathodes 
instead of the coated cathodes and to use narrow-band optical filters -in 

order to see the beam while focusing it. 
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PULSED PLASMA. PROPULSION 

R. G. Jahn and W- von Jaskowsky 
Princeton University, Princeton, New Jersey 


Despite the advanced state of development of ion propulsion devices, and the 
very promising early performance of steady ruagnetop lasmady namic arc thrustors, 
continued interest in the inherently more complex pulsed mode of plasma acceleration 
is sustained by four potentially profitable basic characteristics. Specifically, 
passage of the current through the plasma in intermittent intense pulses may be 
expected (1) to provide higher average thrust density than by steady consumption 
of the same electric power; (2) to involve less heat transfer loss and electrode 
erosion than equal power steady arc devices; (3) to minimize undesirable kinetic 
thermalization within the accelerating plasma; and (4) to generate various unsteady 
electromagnetic field effects which aid in the acceleration and containment of the 
plasma, 

- Practical implementation of these potential advantages requires detailed 
understanding of the participating phenomena. Any pulsed plasma accelerator must 
execute the following controlled sequence of physical processes: (a) the appropriate 

amount -of electrical energy must be accumulated in the external circuit, e.g. in a 
high performance capacitor; (b) a proper distribution of propellant gas must be 
established in the accelerator channel; (c) the circuit must be completed by a 
discharge across the channel electrodes, in the proper geometric location at the 
proper time in the cycle; (d) this discharge must rapidly intensify to a self- 
accelerating current sheet which covers the entire channel cross-section, is 
impermeable to ambient gas ahead of it and its self -magnetic field behind it, and , 
can accelerate over the full channel length stably and with total -ambient gas 
entrainment; ,(e) when the accelerated mass of gas and plasma reaches the end of 
the channel, it must be ejected with a minimum of thermal and electromagnetic loss; 
(f) finally, the entire gas dynamic and electromagnetic system must recover its 
initial situation, prior to initiation of a new cycle. 

Over the past four years our- laboratory has carried out a spectrum of basic 
experiments and related analytical studies encompassing all of the six phases of" 
the pulsed plasma acceleration cycle outlined above. The apparatus has been 
developed primarily from the standpoint of controlled, ' reproducible experimentation 
with maximum diagnostic accessibility, rather than focusing on any particular 
thrustor geometry. At present we have operative five different discharge chambers 
driven by a variety of capacitor banks, pulse-forming networks* and special energy 
storage units. Peak discharge currents up to IQ 6 amp at 10 KV may be attained, 
and pulse lengths from 1/2 to 100 psec are available, diagnostic equipment employed 
includes e.lectric and magnetic probes, voltage dividers, D loops and Rogowski coils, 
streak and Kerr cell cameras, a grating monochrometer and gas laser", high speed 
piezoelectric transducers and 4 mm. microwave bridges and interferometers. Vacuum 

facilities include a 3' dia. x 6' long Plexiglas tank capable of evacuation to 
10"*- 3 torr.\- 3 *Ld0.j 

For purposes of this presentation we will illustrate the program with two of 
the currently active studies: the production of various driving current waveforms 

and their effect on the development of the plasma acceleration process: C'33 ,34) an( j 
the Istribution of current density and magnetic field in exhaust plumes from one 
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type of accelerator. 18 5 20,21 5 29 ,32 ' > With regard to the former we begin by 
noting that while it is -quite conventional to drive pulsed accelerators with 
free-ringing lumped capacitor banks, this is far from an ideal source for this 
purpose. Associated with each reversal of the damped sinusoidal current pattern 
is a secondary "crowbar" breakdown in the discharge chamber which effectively 
uncouples the primary current sheet from the energy source, leaving it to 
dissipate itself over the remainder of its motion down the- channel. (2 , 4,5) 

These secondary breakdowns also badly complicate the experimental diagnosis of 
the discharge-acceleration events by the techniques outlined above. Both 
problems have been eliminated by replacing the lumped capacitor bank with a 
pulse-forming line capable of providing a variety of high current waveforms at 
the operator's discretion. Essentially the device is an L-C ladder network 
formed by two parallel collector plates, along which stock capacitors may be 
attached at arbitrary locations. As suspected, the driving current waveform is 
found to have a profound influence on the effectiveness of the discharge 
acceleration process. Sweeping efficiency and terminal velocity of current 
sheets in a given accelerator can be changed by an order of magnitude by proper 
distribution of the same total electrical energy in the current- pulse. Figure 1 
shows four example current waveforms constructed by this facility and corresponding 
streak photographs of an 8" diameter, linear pinch discharge in ar-gon driven 
by them. 

The importance of the plasma exhaust phase to the over -all performance of 
a pulsed thrustor is evident enough; whatever energy is expended in ejecting the 
plasma slug, or remains trapped within it in various thermal modes, detracts 
from the total thrust efficiency. Although it is difficult to reduce this problem 
to concise elements for systematic study, any information about the development 
of current density and magnetic field patterns in the plume will provide some 
insight into the ejection mechanisms. For this purpose, a special discharge 
chamber with 4" dia. orifice, driven by a 6500 joule pulse line,( 32 ) is allowed 
to eject plasma into a large Plexiglas vacuum tank. 05) Here the development of 
the luminous pattern of the exhaust plume may be determined by Kerr cell 
photography, and the distribution of magnetic fields and current densities 
mapped with magnetic probes and Rogowski loops. Typical patterns are shown in 
FiS* 2. Substantial currents are found to extend far out into the plume, and 
the sense of the local j x B body force is predominantly favorable to continued 
expansion or detachment of the plume. One significant result of such studies is 
the identification of a transient phase of ejection, wherein discrete current 
sheets propagate through the plume, "sweeping" gas along much as they do within 
the discharge chamber, followed by a more steady phasa wherein the current 
patterns stabilize and revert to a "pumping" mode of gas acceleration much like 
that of the steady magnetoplasmadynamic arcs. A second result of interest is a 
demonstrated sensitivity of some aspects of plume development to relatively 
small changes in ambient gas pressure in the exhaust tank, emphasizing the 
importance of test. environment to this phase of the acceleration process. 
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PULSED PLASMA ACCELERATOR 

P, Gloersen, B. Gorowitz, and T. W. Karras 
General Electric Space Sciences Laboratory, King of Prussia, Pa. 


Problem Definition 


This program was initiated for the purpose of developing an efficient 
reliable, lightweight device for ionizing and accelerating gaseous sub- 
stances to velocities the order of 5 x 10 4 meters /second. A coequal part 
of the program is the generation of experimental data which will permit a 
better understanding of the acceleration mechanisms involved in a coaxial 
gun. 

Background 


The motivation for this research is the requirement for thrusters 
with specific impulse characteristics higher than available from chemical 
or nuclear rockets to permit higher payload fractions in interplanetary 
space ^nissions with intermediate characteristic lengths (the order of 
3x10 Km), acceptable payloads for deeper solar system exploration 
(characteristic lengths up to 5x10 Km), and economically feasible solar 
system escape. With respect to competing electric thruster concepts, 
pulsed plasma devices offer the potential of much higher average thrust 
per unit area than electrostatic devices ( py b 2 / 2 y Q vs e G E^/2) and 

much higher peak thrust (with possibly resulting higher coupling efficiency) 
at a given average power level than steady- state MPD arcs. 

Approach 

The program consists of the construction of various accelerator 
configurations based on current understanding of acceleration mechanisms, 
measurement of the gross performance characteristics (thrust, mass flow, 
power input) of such configurations, confirming measurements and detailed 
plasma diagnostics (exhaust stream calorimetry, particle velocity and 
energy probe, magnetic field probe, current probes, particle density probe 
terminal measurements, and spectroscopy) on the most promising, of thes« 
configurations, and the updating of the understanding of the acceleration 
mechanism. (See Figure 1) 
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Results 


In terms of gross performance measurements, the best numbers 
consistently generated to date are an overall efficiency ( r} Q = T^/2m P) 

of 70% or more at' a specific impulse (I Sp = T/mg) of 5000 seconds* 

The power- to-thrust ratio at this point is about 3,6 x 10 4 meters/ second 
(160 KW /pound) and is constant within 5% of this value from 3300 to 5400 
seconds specific impulse. 


The most complete diagnostics to date exist fcr an operating point 
attained earlier, which was f) 0 = 57%, Igp = 4'800 seconds, obtained at a 

mass flow rate of 2*96 x 10" 7 Kg/sec, a capacitance of 144.5 mfd, a 

potential of 950 volts, a pulsing rate of 10 sec , xenon propellant, and 

the A-7D gun configuration (see Figures 2-4). The relations 
* 

2 


7 ] = r? r? 

'o E m 


I = 
sp 


v / g 



2 

cos a 


were used to generate second values for 7 ] 0 and Igp for comparison with 
the results of the gross performance measurements* 7]^ was obtained by 
the use' of a water-cooled calorimeter placed in the exhaust stream during 
repetitive operation of the accelerator under conditions equivalent to those 
cited. An upper limit to ?7 m was obtained by measuring the fraction of 
gaseous propellant available to the discharge by means of a small, fast 
ionization gauge pressure sensor (See Figure 5) used to map the propellant 
density prior to discharge as a function of both time after injection and 

position in the discharge region (See Figures 6-8)„ The factors v, v 2, , 

and cos 2 a were obtained by means of the particle velocity and energy 
probe (See Figures 9-13), taking into account the 7 ] rn of the device. The 

_ 4 

values obtained are 77 • = 72%, T] ~ 90%> v = 4,8 x 10 meters /sec, 

Z ~ xxi 

V / V 2 = 86%, and cos 2 a = 97%. Thus, the agreement between r\ Q and ?? o 

is within 10%, as is the agreement between I S p and I sp which is remarkable 
in view of the estimated cumulative measuring accuracy of 20%. An account- 
ing of current distributions both inside and outside of the accelerator added 
credence to these results, as did the operating background pressure level 
of about 10-5 mm an £ an extended run of 10 hours under equivalent operating 
conditions with no degradation of performance. 


Future Plans 

The mainstream of the next contractual effort will be to adapt the coaxial 
gun to the use of metal vapor propellant and determine the performance of this 
combination. Additional studies using xenon propellant are in progress. 
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Figure 3t, (ay A-7D Accelerator, (b) A-7D Accelerator Central Electrode 
Assembly 







Figure 4. Circuit Diagram of the A7-D Accelerator 
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Gas Density Probe 
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Propellant Mass Fraction vs Value Voltage for Mod A--7D Gun 
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Figure 7>. Fraction of Propellant B.etween the. Accelerator Electrodes 

at Various Discharge Trigger Delay Times 
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NYLON INSULATORS 



Figure 9„ 


Gridded Probe for Measurement of Ion Energies 
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Figure 10. Total Ion Velocity Distribution on Axis, A-7D Gun, 77 = 57*# 
1^ - 4800 sec., 950V, 144.5 \JL fd, .296 mg/sec Xenon 
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Total Ion Velocity Distribution 45 cm. off Gun Axis, A-7D Gun 

77 = 57%, I' = 4860 sec., 950V, 144.5 flid, .296 mg/sec. Xenon 
sp i- 


Figure 'll-*. 
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Figure 1 1 


Fraction of. Accelerated Propellant within a Beam Radius, 2m Downstream from Accelerator 
Muzzle. Mod A-7D, Tj - 57%, I = 4800 sec. , 950V, 144. 5 jUfd, . 296 mg/sec Xenon 
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Figure 13 ^ 


Average Ion Velocity at Various Radial Positions in a Plane Two Meters Downstream From 
the Accelerator Muzzle. Mod A-7D, 7 } ~ 57%, I = 4800 sec., 950V, 144.5 jUfd.> 

.296 mg. /sec Xenon 
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COAXIAL PLASMA GIFT 

A* V. Larson, L. Liebing, A. R. Miller, and R. Dethlefsen 
Space Science Lab., General Dynamics Convair, San Diego, Calif. 

Since 1961, a development program for pulsed plasma guns has been jointly 
sponsored by IIA.SA -Lewis and General Dynamics Convair. The principal Lewis goal 
has been to demonstrate competitive thrust efficiencies in such devices. 'This 
year the number of industry competitors in this field has been reduced to one, 
the decision being based on the efficiency of presently existing accelerators. 

As the Convair program nears completion, it is appropriate at this conference, 
to summarize the effort and results, and to list the related publications. The 
accomplishments of the last year are given in abstracts of recent meetings and 
are reproduced here. Finally, the results wi ll he given of recent experiments 
which were designed specifically to measure the performance of the pulsed plasma 
device of Gloersen, Gorowitz, and Karras. 

The Convair program began with a strong diagnostic effort to examine the 
plasma energy transfer processes. Electromagnetic field prohes were developed 
and applied to measuring the phenomena within the discharge. The techniques of 
Faraday cups, gridded electrostatic analyzers, and other particle analyzers were 
established for the exhausts of pulsed plasma guns. Scaling laws were verified 
for the control of I S p for guns operating with various propellant distributions, 
e.g., uniform fill and slug fill. The applicability of the snow-plow model was 
established. Fundamental limitations, due to the physics of the plasma processes 
in two modes of plasma accelerators were recognized. Important parametric re- 
lationships which affect performance were empirically determined. Three techno- 
logical spin-offs have occurred; an order of magnitude improvement in energy 
storage capacitors; a high speed, low power, compact gas valve; the founding of 
the pulsed MPD arc technology. 

To illustrate our recent work, we include abstracts of two papers given 
at the Fifth Electric Propulsion Conference, March 1966. 

Paper (66- 199 ), " Experiments with a Pulsed Arc Plasma IhrustOr." 

A series of experiments with an improved version of our pulsed arc plasma 
thruster will be described. This thruster is similar to' the MPD arc jets ex- 
cept that it is pulsed. The propellant feed system and the power supply are 
built so that the thruster operates in a steady-state mode for 650 jisec. The 
power level ranges between 10 kw and 1.5 Mw. Hall current, radial current, and 
magnetic field measurements have been made within the thruster discharge. The 
Hall current density is generally an order of magnitude larger than the rad ial 
current density. The jgB r force is about twice larger than the j^Lg force. It 
is estimated that Wt equals 5 Ion the electrons. The thruster current, voltage, 
ion exhaust velocity, and ion density have been measured as functions of the in- 
jected mass flow rate with helium, nitrogen, and argon propellants. The t h r ust 
is not strongly dependent upon the injected Si and, therefore, the apparent effi- 
ciency is very high at low m. It has been found that the cathode erosion becomes 
important at low mass flow rates, but complete erosion studies have’ not been done 

Paper (66-238), "Experiments with a Voltage Switched Pulsed Coaxial Gun." 

An experimental program will be described which investigates the behavior of 
a pulsed coaxial thruster in which slug -model propellant distributions are used. 
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The propellant Is injected radially through ports located near the insulator. 

The gas valve opens in about 100 [isec and a plenum is used to control the 
amount of propellant injected. A triggered voltage switch is -used to delay 
the time between propellant -injection and thruster discharge. Measurements ' 
are made of the neutral -gas distribution prior to the discharge, the current 
distribution, the ion velocity and energy in the exhaust, and 'the calorimetric 
and thrust efficiency. The variables include the firing delay, mass load, 
capacitor voltage, and barrel geometry. The thruster operates in an axially- 
syrsmetric stable manner over an interval in which the gas distribution changes 
from the slug-model approximation to a nearly uniform distribution. Peal calori- 
metric efficiencies of 70 ^ occur for propellant distributions which approximate 
the slug model. The efficiency increases with capacitor voltage and is rather 
insensitive to plenum pressure. 

Since March,, we have concentrated on measuring a, duplicate of the G. E. 
accelerator for which high efficiencies have been reported. Performance has 
been measured as a function of propellant species, port location, direction of- 
propellant injection, firing delay, and barrel geometry. The important result 
is' that it has not been possible in our installation to achieve the efficiency 
reported by General Electric, or to exceed efficiencies achieved by the various 
General Dynamics accelerator designs. 

Publications for MSA Related Work 

-"The use of Magnetic Probes in Plasma Diagnostics,” R. H. Lovberg, Ann. Phys., 

8, 3H (1959). 

"Current Sheet in a .Coaxial Plasma Gun," L. C. Burkhardt -and R. H. Lovberg, Phys* 
Fluids 5, 3^1 (l$62). 

"The Use of a Coaxial Plasma Gun for Space Propulsion," Final Report, Contract 
IAS '5-1139, GP/Astronauties Report Wo. AE 62-0678, dated May 19^2. 

"TJse of Ballistic Pendulums with Pulsed Plasma Accelerators," T. J. Gooding, 

B* Hayworth, and R. H. Lovberg, ARS J;,'Qct. (19&2) p. 1599. 

"Impulsive MED Devices as Space Engines," R. H. Lovberg, Proc. of Third Symp. 
on A‘dv. Propulsion Concepts, held at Cincinnati (Oct. 1962 )', publ. by Gordon and 
Breach, New York, 1963, p. 95* 

"Physical Processes in a Coaxial Plasma Gun," T. J. Gooding, R* R. Hayworth, 

R. H. Lovberg, AIM. Electrical Propulsion Conference, Paper 63-004, Colorado 
Springs, March ( 1962 )# 

"Instabilities in a Coaxial Plasma Gun," T. J. Gooding, B, Hayworth, and R. H. 
Lovberg/ AIAA J. 1, 1289 ( 1963 ). 

"Development of a Coaxial Plasma Gun for Space Propulsion, " Final Report, 

Contract MS 3-2501, GD/Astronautics Report Wo. GDA £>3-0454, dated May 1963. 

"inference of Plasma Parameters from 'Measurement of E and B Fields in a Coaxial 
Accelerator," R. Lovberg, Phys. Fluids 7, Part 2-Supplement^ S57 (l964) . 
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"development of a Coaxial Plasma Gun for Space Propulsion, " Final Report, 
Contract -NAS 3-2594, MSA. Report, No.. CR-54l49, GD/Astronautics- Report No. 

GBA DBE 64-051, dated June 1964.. 

"An Energy Inventory in a Coaxial Plasma Accelerator Driven "by a Pulse Line 
Energy Source," A. V. Larson, et al., AIAA J.> -3,- 977 (1965) • 

"Energy Loss in Pulsed Coaxial Plasma Guns," D. E. 1. F. Ashby, AIAA J. 3,- 

1045 (1965). 

"Exhaust Measurements on the Plasma from a Pulsed Coaxial Gun," D. E. T. F .. 
Ashby, 5. J. Gooding, B. R. Hayworth, and A. V. Larson, AIAA J. 3>, 1140 ( 1965 ); 
also AIAA Reprint 64-704 (1964)'. 

- t 

"Development of .a 'Coaxial Plasma Gun for Space Propulsion," Final Report, 
Contract MS' 3-5759, MSA Report No. CR-54245, dated April 1965 . 

"Energy Storage Capacitors for Pulsed Plasma Thrustors, " B. R. Hayworth, 

A. R. Miller./ and C. U. White, AIAA Paper No. 65 - 337 , San Francisco, July 1965- 

/ 

"Quasi-Steady-State Pulsed -Plasma Thrustors," D. E. T. F. Ashby,- T. J. Gooding, 

and A. V. Larson, AIAA Paper No. 65 - 338 , San Francisco, July 19§5; to. be 

published (AIAA J.). 

\ 

"Note on the Use of Mass Analyzers in Plasma Physics," H. Fleischmann, 

D. E. 1. F. Ashby, A. V. Larson, (to be published in Nuclear Fusion); also 
available as Report GA-6219, General Atomic, San Diego, Calif. 

"Generation of Plasm Streams with Pulsed Arc Discharges," T. J. Gooding, 

D. E. T. F. Ashby, and L. Liebing, General Dynamics Convair Report No. ERR-AN 
739, May 1965 . 

"Experiments with a Pulsed Arc Plasma Thrustor," L. Liebing, A. V. Larson, 

T. J. Gooding, AIAA Paper No. 66-199, Fifth Electric Propulsion Conference, 

San Diego, March 1966. 

"Experiments with .a Voltage .Switched Pulsed Coaxial Gun," A. V. Larson, 

L. Liebing, A. R, Miller, AIAA Paper No. 66 - 238 , Fifth Electric Propulsion 
Conference, San Diego, March 1966 . 
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DIAGNOSTICS OF ACCELERATING PLASMA 

Drs. C. C. Chang, T; N. Lie. and A.W.. Ali 
The Catholic University of America 

ABSTRACT 


In pulsed plasma accelerator, the process-by which the neutral gas is ionized,' 
heated and acquires momentum is not fully understood. A simple model in which a 
shock wave, with no losses, is accelerated by magnetic piston is certainly inadequate-. 

In. order to understand these processes, the transient plasma state must be known, i.e., 
electron densities/composition, temperatures are required for comparison with theoret- 
ical models. 

1 

' In our project, it is intended: a) to apply spectroscopic techniques to diagnose 
the plasma produced in some typical plasma accelerator and compare the results with 
those obtained by other investigators; b) to compare existing theory for plasma acceler- 
ation with the measurement; c) to account for the energy deposited in the plasma in order 
to give further information on losses. 

As a first step of the project, accelerated plasma produced by a coaxial plasma 
gun which was operated at static ambient pressure (helium) was studied. ‘Electron tem- 
perature; electron density, and other particle densities are determined at a location 
which is 1 cm outside the gun end. Furthermore, the electron temperature and the 
density were obtained as a function of radius of outer electrode using Abel inversion 
procedure with a series of measurements. 

The same technique has been extended into the gun barrel to diagnose accelerating 
current layer. The electron densities and temperatures inside the gun were relatively high 
and were found to be Te~ il eV and Ne — 7 X 10^7 C m - 3 at the location which is 
7.5 cm from the insulator. 
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‘MED ARC 'JET TRRUSTOR RESEARCH 

Contract HAS 3-2593 
AVCO Corporation 
Space Systems Division 
Wilmington, Massachusetts 

Principal - Dr. R. R. John and Ur. S. Bennett 

_ Associates - Mr. .G. Enos and' Dr. A-. Tuc hma n 

The overall purpose of the subject research and development 
program is to establish whether the laboratory MPD arc jet, which appears to "be 

an efficient gas accelerator, can be converted into a reliable low thrust 'space 
propulsion engine. 


Investigators : 


Objective: 


| teckgrou n d ^ ^ _ A typical MPD arc jet configuration is shown in figure 1. 

The basic unit is axially symmetric, consisting of a central cathode surrounded 
by a co-axial anode, and can be operated both with and without a magnetic field 
The working fluid is accelerated both by aerodynamic and magnetoplasmadynamic 
forces as it passes through and around the discharge. Although there is con- 
siderable controversy regarding the exact nature of the gas acceleration 
process, the basic thrust producing mechanisms characteristic of the MPD arc let 
associated with: (l) Aerodynamic pressure forces (thermal abceleration) : 

( ) Self- induced magnetic field forces; and (3) External magnetic field forces. 


Progress and Results ; Five basic thrust producing mechanisms- have been pro- 
posed as being characteristic of plasma generators of the type shown in figure 
1. These are: (1) Aerodynamic pressure forces; (2) Magnetic pumping; (3) 

- Ma f n ^ c blowin S^ (*0 Aerodynamic swirl introduced by electromagnet forces; 
aa ;J2 Hal1 cm>rent acceleration. Aerodynamic pressure forces are associated 
with thermai expansion of a gas within confining solid walls; magnetic pumping 
Jesuits' from the interaction of the axial discharge current density, j = , with 

d / Zim S hal m S netic field > B e, i-e. T pump = l/2 I?; magnetic 
lowing results from the interaction of the radial component of the current 

r2 n /W)!\ 3 F* W1 ^ h t* 16 self -i»duced azimuthal magnetic field, B a , i.e. Thi ow = 

I (1/4- +y.n r a /rc)., where r a and r c are respectively the radial extent of 
he current distribution at the anode and cathode. These three thrust mechan- 
isms have been clearly identified and their existence is not dependent on the 
presence of a magnetic field?- The remaining mechanisms are associated with 
an external field. The magnetic swirl is associated with the bulk rotation 
of the gas and results from the interaction of j r with B a . Finally, the Hall 
current term is associated with the interaction of j e with B r . Ho simple ' 
analytical expressions have been derived for the external field terms except 
empirically it has been established that - T " 

f (i,b) . total aero Tse lf " ^external = 


As was originally pointed out by Fahleson for the case of a plasma rotating 
a magnetic field , and by Cann for the MPD arc jet3, the available gas flow, 
k^and the .gas flow, mj, accelerated and ionized-by 'the discharge are not, a 
priori, equal. The concept of a characteristic discharge mass flow has been 
strikingly, verified experimentally at this laboratory by the observation that 
the engine will operate and produce thrust with no mass flow passing through 
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the propellant feed line 11 -. - -By. analogy with conventional arcs' Caim3 proposed 
that the MPD discharge would operate at a power level such that the Input 
power is .minimized; Fahleson 2 / on the other hand, following a suggestion of 
Alfven, postulated the discharge would operate, .at a mass flow such that the 
kinetic and ionization energies of the particles in the discharge mass flow 
are equal. Assuming that the input power- is primarily absorbed in kinetic 
and ionization energy the two postulates are 'identical. From energy consider- 
ations the flow velocity 'at minimum input power is 'given by Ucrit = -\/2eYi/ 
maj the critical mass flow becomes md,min * t/ \/2eYi/m a ; and the minimum power 
is IV) m i n = T^eV i/t%. The empirically observed linearity between Y and B 
and T and IB, shown in figures 2 and 3, is compatible with these minimum 
energy considerations. 

-The preceding discussion has been based on a macroscopic ■ description of accel- 
erator performance. It has been proposed that the external magnetic field 
forms a magnetic nozzle in which azimuthal kinetic energy resulting from the 
interaction of B a and j r is !t spun out” into axially directed kinetic energy; 
however, the exact nature of acceleration and confinement process is hot' at 
all clear. Local current measurements in the MPD axe jet discharge obtained 
with Hall effect sensors, suggest that the basic acceleration process occurs 
within two or three diameters of the engine 'face. The local cathode tip 
pressure is observed to decrease with increase in magnetic field suggesting 
magnetic confinement effects; and, finally, local velocity measurements with 
a Faraday probe indicate strong radial velocity dependence. The local diag- 
nostic measurements have, however, not yet been carried out in sufficient 
detail to distinguish between different possible acceleration mechanisms. 

Plans; * . The future direction of the research program will he to 

continue to attempt to clarify the nature of the acceleration process by means 
of local diagnostic measurements. The Doppler shift technique will be used to 
establish both axial and' azimuthal ion velocities within the discharge region. 
Detection of significant azimuthal kinetic energy would imply that the princi- 
ple thrust production mechanism is indeed due to j r B g interaction. From the 
engine development viewpoint effort will be directed to the design and con- 
struction of a -radiation-cooled MPD arejet engine for the power range 5 to 50 
kilowatts . 


References , ■ 

1. John, K. R., Bennett,. S., ana Connors, J. F., "Experimental Performance of 
a High Specific Impulse. Arejet Engine", Astronautica- Acta 11, 9T (1965). 

2. -Fahleson, U. Y. , "Experiments with Plasma Moving Through Heutral Gas", The 
Physics of Fluids, Vol. 4, Ho. 1, pg 123-12T (1961). 

3. Gann, G. L. et al, "High Specific Impulse Thermal Arejet Thrustor Tech- 
nology", AFAPL-TB-65-48 Part I, AF Aero Propulsion Laboratory (1965). 

4. John, R. R. , and Bennett, S., "Recent Advances, in Electrothermal and Hybrid 
Electrothermal-Electromagnetic Propulsion, AFOSR Fourth .Symposium on Advanced 
Propulsion Concepts, Palo Alto, Calif., (1965). 


Additional -Papers 

Powers, W. E., Measurements of the Current Density Distribution in the Exhaust' 
of ah MPD Arejet", AIAA 3rd Aerospace Sciences Meeting, AXAA Paper Ho. 66-116, 
January, 1966. ' 


Bennett, S., Enos, G.., -John, R. R., and Tuchman, A., "Cesium Fueled., MPD Arejet 
Engine Performance", 2nd Annual AIAA Meeting, 1 AIAA Paper Ho. 65-296, July 1965. 
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Figure 1 SCHEMATIC DIAGRAM OF THE HYBRID 
(ELECTROTHERMAL-ELECTROMAGNETIC) ARC JET 





FIGURE 2 ARCV&LTAGE VERSUS APPLIED AXIAL 
MAGNETIC FIELD FOR AN HYDROGEN- 
FUELED MPD ARCJET 





FIGURE 3 ENGINE THRUST VERSUS PRODUCT OF 
ARC CURRENT AND APPLIED AXIAL MAG- 
NETIC FIELD FORA HYDROGEN -FUELED 


MPD ARCJET 
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9th NASA’ Inlercenler and Contractors' Conference on Plasma Physics 

Washington, D! C. , May 24-27, 1966 

PERMANENT MAGNETS' FOR MPD ARC THRUSTERS 
(Contract NAS 3-7923) 

A.. C. Eckert and'D. B. Miller 
General Electric Space Sciences Laboratory 
King of Prussia, Pennsylvania 

ABSTRACT - 

A cylindrical Alnico 8 permanent magnet. assembly is being .tested, 
in order, to evaluate the magnet requirements to satisfy MPD thruster 
-field needs. A multi's actional design has been employed so that the effect 
of magnet length and .area oh gap field can be evaluated. Gap field 
strengths of over 2000 gauss have been generated. The measured 
dependence of gap field on magnet geometry will be discussed. 
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Washington, D. C. , May 24-27, 1966 

ELECTRON CYCLOTRON RESONANCE 
PLASMA THRUSTER DEVELOPMENT 
(Contract NAS 3-8903) 

D. B. Miller, A. C. Eckert and G. S. Cook 
General Electric Space Sciences Laboratory 
King 'Of Prussia, Pennsylvania 

ABSTRACT 

The short version of the X-band ECR accelerator has been operated 

on a thrust stand at c-w r-f power levels up to 2000 watts. Although thrust 
2 

efficiencies (T /2mP) of over . 30 have been measured, a good deal of 
randomness has yet to be understood. A retarding potential ion energy 
analyser has been employed simultaneously with the thrust stand to obtain 
an independent measure of specific impulses and velocity distribution. Ion 
probe specific impulses are significantly lower than the corresponding 
thrust-based I , possibly as a result of the relatively important charge 

Sp v 

exchange process occurring with these experiments. Measured velocity 

— 2 2 

distribution efficiencies ((v) /v ) are generally greater than . 80. 
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HALL CURRENT ACCELERATOR 

G. L. Cann, P. F. Jacobs 
Electro-Optical Systems, Inc., 

A. Subsidiary of Xerox Corporation 

ABSTRACT 

A program to investigate experimentally and analytically axisymmet- , 
ric Hall current accelerators has been completed. The results are reported 
in Ref. 1. The purpose of the investigations was to determine the perfor- 
mance potential of such a device as an electric space propulsion engine. 
Parametric studies of the accelerator were made with several propel lants, 
emphasis being placed upon the tests using hydrogen and sodium. The 
accelerators used in these tests are shown schematically in Figs, la and 
lb. The overall performance measured in these devices is presented in 
Figs. 2a and 2b. Some diagnostic measurements on the engine and exhaust 
beam were made and are reported, the most important being a determination 
of the current path in the exhaust beam of the hydrogen accelerator. The 
results of one such measurement is shown in Fig. 3. Mechanisms to explain 
thrust production and other engine characteristics are postulated and 
discussed in detail. The concept of an '’effective" mass flow rate, which 
is determined by a minimum potential hypothesis, is developed to the 
point where a complete analysis of a simplified model of the accelerator 
can be made. Such an analysis is completed and the results are compared, 
where possible, with experimental results. 

Adequate detailed knowledge of accelerator mechanisms to use in an 
engine optimization program was not generated by the above program nor is 
it available elsewhere. A detailed diagnostic program, coupled with an 
analytic investigation, has been undertaken to fulfill this need in our 
present contract work with NASA. 

The fundamental goals of this study are to investigate the opera- 
tion of an alkali metal accelerator with background pressures less than 
10“^ torr, input power 5 -* 50 kw over the specific impulse range of 1000- 
6000 sec. A parametric investigation of overall performance (thrust, 
thrust efficiency and thermal efficiency) with variations in input power, 
mass flow rate and magnetic field is to be carried out. The following 
diagnostic measurements are to be made: 

1) Axial Ion Velocity 

2) Rotational Ion Velocity 

3) Thrust 

4) Torque 

5) Magnetic Field Distribution 

6) Current Density Distribution 

7) Total Beam Power 

8) Local Mass Flux 

9) Local Energy Flux 
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The axial ion. velocity is being measured by both Doppler spectros- 
copy and tracer techniques. The ion rotational velocity is being 
determined solely by spectroscopy. The accelerator thrust is measured 
with a reaction balance, and the angular momentum in the beam will be 
measured with a torque target. The spatial distribution of the three 
components of the magnetic field will be measured with a Hall-effect 
magnetic probe composed of three mutually orthagonal sensing elements. 

The distribution of current density within the plasma is determined from 
these measurements through graphical differentiation and the use of 
Maxwell's equations. The total beam power will be measured with a 
calorimeter composed of 16 separate, water cooled, concentric, ring 
elements. The use of separate elements is important not only from the 
standpoint of a more detailed distribution throughout the plasma, but 
also to avoid shorting the radial electric fields which are probably 
present. The local energy and mass flux distributions will be measured 
with two new probes. The former consists of an open ended, water cooled, 
calorimeter and the latter depends upon the time varying resistance of a 
deposited layer of alkali-metal (viz. lithium) between two collecting 
elements on an insulating substrate. The time rate of change of the 
current is then related to the instantaneous local value of the alkali 
metal mass flux. 

Preliminary results indicate rotational ion velocities of about 
4000 meters/ second at a radial position 1 cm from the jet centerline. A 
number of spectral lines associated with transitions between the upper 
excited states of the lithium ion have definitely been observed. Since 
these states have energies in excess of 50 electron-volts, relative to 
the atomic ground state, this indicates reasonably high plasma electron 
temperatures. The preliminary tests with the mass deposition probe 
indicate the utility of this device as a plasma diagnostic tool. Further 
calibration and elimination of effects related to variations in the local 
plasma potential are still necessary. . 


Reference ; 


1 . 


G. L. Cann, et al., "Hall Current Accelerator" NASA Report CR- 
54705, 4 February 1966 (Contract NAS 3-5909, NASA-LRC) 
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FIG. la HALL CURRENT ACCELERATOR WITHOUT TUNGSTEN ANODE INSERT 

(h 2 -i) 




112 














EFFICIENCY, rf 0 






16 20 24 28 


cm 


JETS AS A FUNCTION OF AXIAL POSITION 



PRECEDING PAGE BLANK NOT FILMED 1 


117 


PLASMA ASPECTS OF SPACE PHYSICS 

A. G. Opp /Physics and Astronomy Programs 
Office of Space Science and Applications, NASA 


One of the objectives of the Physics and Astronomy Programs is to 
understand the physical processes occurring in space. Space is one 
of the most nearly perfect low density collisionless plasma laboratories 
available. Thus, many of the phenomena which are of interest to space 
physicists are also of great interest to plasma physicists, and vice 
versa. 

Space plasmas exist from the lower boundary of the ionosphere to the 
outer limit of the Universe. Before the advent of satellites the 
terrestrial ionosphere was probably the best understood space plasma. 
None the less, ionospheric satellites enabled observations to be made 
above the range of ground based ionosondes and increased immeasurably 
the knowledge of the topside of the ionosphere. They also presented 
plasma physicists with observations of a new phenomenon, the so-called 
Alouette resonances. 

Beyond the ionosphere one finds the trapped radiation belts. The 
mechanisms of trapping, particle acceleration and containment in the 
earth s magnetic field are still poorly understood, as is the magne- 
tospheric tail and its relation to the radiation belts and auroral 
zones . 

The earth's magnetic field and the radiation belts are terminated at 
approximately 10 earth radii in the subsolar direction by the solar 
wind. ~ The transfer of energy from the solar wind to the magnetosphere 
is one of the most challenging problems of contemporary space plasmas. 
Extensive satellite and theoretical investigations of the interaction 
of the solar wind with the earth! s field are being conducted in an 
attempt to understand the physical processes associated with the energy 
transfer. 


The nature and origin of the solar wind is also being actively investi- 
gated. The existence of the solar wind has been firmly established. 

It remains, however, to define its detailed spectral characteristics 
and temporal variations, particularly as these can be related to events, 
on the sun or in near interplanetary space. 

The sun has been studied for centuries. Yet there, still exist no fully 
satisfactory theories for many solar phenomena. Undoubtedly, the ultimate 
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explanations will have their foundations in , plasma physics, and 
.conversely, solar physicists will undoubtedly present plasma physicists 
with new phenomena for explanation. 

The radial extent of the solar wind is a matter of speculation. 
Estimates range from 2 or 3 to over 100 astronomical units. The inter- 
action of the solar wind with the galactic cosmic ray gas and galactic 
magnetic field are matters of considerable interest to cosmologists as 
well" as the origin and propagation of the cosmic ray gas in the Galaxy 
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' EXPERIMENTAL AND THEORETICAL INVESTIGATION ON 
SELECTED ASPECTS OF PLASMA TURBULENCE 

L Bo" Bernstein, Yale University; New "Haven, Conn , 

The theoretical work accomplished thus.- far has included 
studies of the quasi-linear theory of plasma waves appropriate for . 
an understanding of plasma turbulence, and an analysis of the sta- 
bility of the Chapman -Ferraro model of the magnetopause against 
gusting in the solar wind. These latter studies are being extended 
using an ideal magnetohydro dynamic model which includes the bow 
shock.. Also under way are studies of the interaction of electro- 
magnetic waves and the proper oscillations of the plasma in systems . 
exhibiting sheaths. 

Experimental studies already well under way are in the 
areas of (a) light scattering from plasma fluctuations, (b) study 
of cyclotron resonance plasma echoes, and (c) study of finite geo- 
metry plasma resonances; In (a) both plasma production and laser 
stabilization are essentially solved; in (b) collisionless echoes 
due to relativistic mass shift have been observed; in (c) resonance 
spectra both with and without streaming are undergoing study and 
analysis . 
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INVEST IGATIOH OF PLASMA RESONANCE PHENOMENA 
S. J. Tetehbaum 

Lockheed Research Laboratories, Palo Alto, California 


This investigation was undertaken to obtain. an understanding of electron 
resonance phenomena in laboratory and ionospheric plasmas and to evaluate the 
applicability of such phenomena to the measurement of charged particle den- 
sities, temperatures and magnetic field strengths in planetary atmospheres 
and interplanetary space. Experiment's were performed on RF-excited, elec- 
trodeless, pulsed and steady state, low-pressure Helium discharges of 
relatively simple geometry which simulate a plasma slab in a uniform magnetic 
field. The wave frequency, w, the magnetic field intensity, B 0 = mw c /e, the 
gas pressure, p, the electron density, electron temperature and the geometry 
were chosen to provide an approximate correspondence between the laboratory 
plasma and the plasma environment of the Alouette satellite. The major 
resonances studied were the electron cyclotron harmonic resonances which 
occur near w = Nw where N is a positive integer. A schematic of the experi- 
mental arrangement is shown in Figure 1. The plasma was created by a 22 MHz 
pulsed or CW transmitter capacitively coupled to a rectangular Pyrex vessel 
having a 16,5 x 16.5 cm cross-section and an inside thickness d of 9.7 cm in 
the direction of propagation. Its walls were ground to a thickness of a 
half -wavelength at 9-15 GHz, the frequency of the incident wave P^ n . The 
vessel was located in a DC magnetic field which was uniform to better t han 
0 . 5 $ over the plasma volume. The measured magnetic field positions of the 
resonances are estimated to be accurate to 0.1^, A 36 . 0 GHz interferometer 
measured the average electron density over a 2 cm diameter cylindrical 
region through the plasma coaxial with the interferometer horns. Measure- 
ments were taken of P , the incoherent noise emission, of (P w )m and (P w )g, 
the coherent transmitted and reflected fundamental waves and ox (Pgoi)^ and 
' E 2 w)r> corresponding second harmonic waves. The CW microwaves propa- 
gated in free space on either side of the plasma container in a direction 
perpendicular to the static magnetic field. The waves were launched and/or 
detected hy microwave horns located a few wavelengths from the plasma. The 
horns could be rotated to propagate or receive either the extraordinary wave 
(%. 15-fBo) ° r tlie ordinary wave (Eq jn |j B Q ) . The intensity of a given wave 
at a fixed frequency of 9*15 GHz or 18.3 GHz was recorded as the magnetic 
field was slowly swept in intensity and precise measurements were made of 
the positions, amplitudes and line shapes of the resonances under varying ' 
conditions of the plasma parameters. Figure 2 shows the resonance amplitudes 
of some of the waves as a function of harmonic number 13, where smooth curves 
have been drawn through the measured points. The incoherent wave resonances 
reach maximum amplitude near N = 4 or 5 and thereafter gradually decrease 
with increasing 13, the ordinary wave resonances being 1 or 2 dB less than 
the extraordinary wave resonances. The latter were observed up to N = 46. 

The coherent wave resonances were only observed when the incident microwave 
electric field was perpendicular to the magnetic field. Fundamental wave 
resonances were observed up to N = "J > their amplitudes decreasing with 
Increasing N . They exhibited an increase in power relative to the back- 
ground power in transmission, and ,a decrease in power relative to the 
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background power in reflection. Second harmonic wave resonances were observ- 
ed up to W 7, They exhibited an increase, in power relative to the back- 
ground power in, both transmission and reflection. The amplitudes of the 
(P 21 o)t resonances were about 3 dB greater than those of (3?g u )p; and in some 
cases were 1 6 cLB above the second harmonic wave background . The rapid 
decrease in resonance amplitude with increasing U may be related to the 
rapid decrease in the generated second harmonic in this magnetic field 
region. Resonances were also observed in the sum frequency wave., „ P S' due to 
two incident waves of frequency Wj/W = 8.8 GHz and « 9.2 GHz. These 

resonances occurred close to w^. - w 2 " ^ w c-» ^ = 2, 3 , 4 and 5 • 

Attempts to observe resonances in P 2w near 2w = Hu , H ~ 3*5*7* and in Pg 
near = Bs) c , H = 3,4,5. . • were unsuccessful. The coherent wave and 

incoherent wave resonances exhibit different behaviour not only in their 
amplitude dependence but also in their positions. Figure 3 shows how the 
emission resonance positions vary with N. The positions remained constant 
for all harmonic numbers, for pressures from 0.02 to 0.3 Torr and for values 
of oup/u) at the center of the discharge from 0.8 to 2.6, where w p is the 
electron plasma frequency. The resonance positions of the two waves 
differed by 0.3$, those of the extraordinary wave occurring at a higher 
magnetic field than those of the ordinary wave. In contrast, the coherent 
wave resonances were complicated functions of If, p and m /to. Their positions 
varied as much as 4$ with variations in these parameters, but always occurred 
on the .low frequency side of the exact harmonic . The coherent and 
incoherent wave resonance amplitudes also differ in their pressure depend- 
ence. The former show definite collisional damping effects, the amplitudes 
decreasing with increasing pressure, while the latter exhibit a maximum 
amplitude at a pressure whose value depends upon m /to, the amplitude 
decreasing for higher and lower pressures. The incoherent wave resonances 
had symmetrical but non-Lorentzian line shapes. The Q’s of the extra- 
ordinary wave resonances' were about 40 and were about one third larger than 
the Q/s of the ordinary wave resonances. The coherent wave resonances had 
asymmetrical line shapes, their asymmetry decreasing with increasing 1,. and 
Q’s of about 30. The coherent wave resonance lines were broader on their 
high magnetic field side and occasionally, for 3f = 2 and Iff = 3* showed the 
presence of additional small resonances. Experiments in which the center of 
the plasma was moved relative to the center of the magnetic field have 
shown that the cyclotron harmonic emission is a volume effect father than a 
sheath effect and that the radiated ordinary and extraordinary wave 
resonances originated in different regions of the plasma. 
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A-C TRAVELLING WAVE-PLASMA STREAM INTERACTIONS 

Martin Lessen 

The University of Rochester 
Rochester, New York 

NsG-350-63 

The motivation of this problem was to investigate the 
interaction of a travelling magnetic wave with a plasma. 

The configuration selected for this research was a pair of 
coaxial cylinders with an inner and outer shell of low re- 
luctance magnetic material. The configuration was chosen 
to provide an electrodeless machine with a nearly radial 
magnetic field. The concept of this machine was first 
presented in reference (1) in which a single fluid model 
was used for analysis. Subsequently, this type of confi- 
guration was discussed by Carter and Laubenstein (2) and 
more recently by Penfold et al (3) and Schwirian (4). 

The approach was initially in two directions. The 
first was the analysis of a three fluid one dimensional 
model of the device in which the principal mechanism for 
plasma neutral interaction was assumed to b’e charge ex- 
change. The second was experimental which involved con- 
struction of the coaxial cylindrical device including the 
windings low reluctance, magnetic material shells, gas 
metering system and some diagnostic elements. This latter 
area indicated that a great deal had to be done before sig- 
nificant physical results could be obtained as to the in- 
ternal interaction mechanisms. 

Initially several sets of windings were made and tuned 
in an attempt to obtain the strongest radial field with as 
little axial component as possible. The final design con- 
sisted of four 9 turn coils each approximately 9 inches 
long. For this configuration the radial field strength 
varied from about 100 gauss to 160 gauss along the tube - 
having maximums between the coils. At the same time, the 
axial component was about 20 to 25 gauss representing an 
angularity of the field midway in the annulus of some 
significance. The phase angle of the field fluctuated 
over a 12.5° range along the tube from an assumed linear 
360° total phase shift. Due to the field angularity the 
reaction with the plasma would drive many of the particles 
to the walls before the exit were reached, however, there 
should be some measurable reaction. The most direct measure- 
ment was of the pressure distribution along the tube for a 
given flow rate. By having the field accelerate the plasma. 
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a pressure drop from 300yiHg to IlOyuHg was recorded over 
3/4 of the tube length; connecting the field in reverse 
necessitated an inlet pressure of approximately 800/tHg 
to realize the lOOy^Hg at the same point. Considerable 
changes in the electrical power input were noted during 
these interactions. Some estimate of the power transferred 
to the plasma was desired and coils wound in the annulus and 
connected to non-inductive resistances were used for this 
purpose. Induced power, in -the coils amounted to approximately 
6 KW. At this time the power to the coils totaled 10.2 KW. 

In the discharge mode tjie power dropped to 5 KW to the coils 
due to considerably more detuning. 

Work is continuing with the construction of a new 
machine with conical geometry which should alleviate some 
of the problems of field angularity and the small annulus. 

Also more room will exist for diagnostics than could be 
affected in the current machine. It is anticipated that 
the question of the magnitude of elastic scattering of 
particles and radial electric field effects can be looked 
at as well as the continued investigation of overall in- 
teraction. 
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